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Abstract 



We introduce a mean-reverting SDE whose solution is naturally defined on the space of correlation 
matrices. This SDE can be seen as an extension of the well-known Wright-Fisher diffusion. We provide 
conditions that ensure weak and strong uniqueness of the SDE, and describe its ergodic limit. We also 
shed light on a useful connection with Wishart processes that makes understand how we get the full 
SDE. Then, we focus on the simulation of this diffusion and present discretization schemes that achieve a 
second-order weak convergence. Last, we explain how these correlation processes could be used to model 
the dependence between financial assets. 
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Introduction 

The scope of this paper is to introduce an SDE that is well defined on the set of correlation matrices. Our 
main motivation comes from an application to finance, where the correlation is commonly used to describe 
the dependence between assets. More precisely, a diffusion on correlation matrices can be used to model 
the instantaneous correlation between the log-prices of different stocks. Thus, it is also very important for 
practical purpose to be able to sample paths of this SDE in order to compute expectations (for prices or 
greeks). This is why an entire part of this paper is devoted to get an efficient simulation scheme. More 
generally, processes on correlation matrices can naturally be used to model the dynamics of the dependence 
between some quantities and can be applied to a much wider range of applications. In this paper, we mainly 
focus on the definition, the mathematical properties and the sampling of this SDE. However, we discuss in 
Section |4] a possible implementation of these processes in finance to model a basket of risky assets. 

There are works on particular Stochastic Differential Equations that are defined on positive semidefinite 
matrices such as Wishart processes (Bru ) or their Affine extensions (Cuchiero et al. [8 ) . On the contrary, 
there is to the best of our knowledge very few literature dedicated to some stochastic differential equations 
that are valued on correlation matrices. Of course, general results are known for stochastic differential 
equations on manifolds. However, no particular SDE defined on correlation matrices has been studied in 
detail. In dimension d = 2, correlation matrices are naturally described by a single real p E [—1, 1]. The 
probably most famous SDE on [—1, 1] is the following Wright-Fisher diffusion: 
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where k > 0, p E [—1,1], cr > 0, and {Bt)t>o is a real Brownian motion. Here, we make a slight abuse 
of language. Strictly speaking, Wright-Fisher diffusions are defined on [0, 1] and this is in fact the process 
> 0) that is a Wright-Fisher one. They have originally been used to model gene frequencies (see 
Karlin and Taylor 20 ). The marginal law of Xt is known explicitly with its moments, and its density can 
be written as an expansion with respect to the Jacobi orthogonal polynomial basis (see Mazet [13]). This 
is why the process {Xt,t > 0) is sometimes also called Jacobi process in the literature. In higher dimension 
(d > 3), no similar SDE has been yet considered. To get processes on correlation matrices, it is instead used 
parametrization of subsets of correlation matrices. For example, one can consider Xt defined by {Xt)ij = pt 
for 1 < i j < d, where pt is a Wright-Fisher diffusion on [—l/{d— 1), 1]. More sophisticated examples 
can be found in [21] . The main purpose of this paper is to propose a natural extension of the Wright-Fisher 
process ^ that is defined on the whole set of correlation matrices. 

Let us now introduce the process. We first advise the reader to have a look at our notations for matrices 
located at the end if this introduction, even though they are rather standard. We consider {Wt,t > 0), a 
d-hy-d square matrix process whose elements are independent real standard Brownian motions, and focus 
on the following SDE on the correlation matrices £^(1^): 

Xt=x + [ (At(c - Xs) + (c - Xs))i) ds + V an / (V^^ - ^se^X.dVF.eS + e^W]^ ^X., - X^e^X^) , 
Jo Jo ^ ' 

where x,c € ^d(M-) and k = diag^Ki, . . . , Kgi) and a = diag(ai, . . . ,ad) are nonnegative diagonal matrices 
such that 

KC + CK-{d-2)a'^ eSj{R) or d^2. (3) 

Under these assumptions, we will show in Section [2] that this SDE has a unique weak solution which is 
well-defined on correlation matrices, i.e. Vi > 0,Xt £ £rf(K.). We will also show that strong uniqueness holds 
if we assume moreover that x S £^(R) and 

Kc + CK-da^ eSj{R). (4) 

Looking at the diagonal coefficients, conditions ([3]) and (jU) imply respectively Ki > {d ~ 2)af/2 and > 
daf /2. This heuristically means that the speed of the mean- reversion has to be high enough with respect to 
the noise in order to stay in £d(R). Throughout the paper, we will denote MRC'd{x,K,c,a) the law of the 
process {Xt)t>o and MRCd{x,K,c,a;t) the law of Xt- Here, MRC stands for Mean-Reverting Correlation 
process. When using these notations, we implicitly assume that ([3]) holds. 

In dimension d = 2, the only non trivial component is {Xt)i^2- We can show easily that there is a real 
Brownian motion {Bt,t > 0) such that 

diXt)i,2 - (ki + K2)(ci,2 - iXt)i.2)dt + ^al+al^l-{Xt)l^^dBt. 

Thus, the process ([2]) is simply a Wright-Fisher diffusion. Our parametrization is however redundant in 
dimension 2, and we can assume without loss of generality that ki = k.2 and ai — a2. Then, the condition 
nc + CK e 5^^(11^) is always satisfied, while assumption ([4]) is the condition that ensures > 0, (^4)1^2 G 
(—1,1). In larger dimensions d > 3, we can also show that each non-diagonal element of ([2]) follows a 
Wright-Fisher diffusion ([Ij. 

The paper is structured as follows. In the first section, we present first properties of Mean-Reverting 
Correlation processes. We calculate the infinitesimal generator and give explicitly their moments. In partic- 
ular, this enables us to describe the crgodic limit. We also present a connection with Wishart processes that 
clarifies how we get the SDE ([2]). It is also useful later in the paper to construct discretization schemes. Last, 
we show a link between some MRC processes and the multi-allele Wright-Fisher model. Then, Section [5] is 
devoted to the study of the weak existence and strong uniqueness of the SDE ^ . We discuss the extension 
of these results to time and space dependent coefficients k, c, a. Also, we exhibit a change of probability 
that preserves the family of MRC processes. The third section is devoted to obtain discretization schemes 
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for This is a crucial issue if one wants to use MRC processes effectively. To do so, we use a remark- 
able splitting of the infinitesimal generator as well as standard composition technique. Thus, we construct 
discretization schemes with a weak error of order 2. This can be done either by reusing the second order 
schemes for Wishart processes obtained in [5] or by an ad- hoc splitting (see Appendix [D|. All these schemes 
are tested numerically and compared with a (corrected) Euler-Maruyama scheme. In the last section, we 
explain how Mean-Reverting Correlation processes and possible extensions could be used for the modeling 
of d risky assets. First, we recall the main advantages of modeling the instantaneous correlation instead of 
the instantaneous covariance of the assets. Then, we discuss our model and its relevance on index option 
market data. 



Notations for real matrices : 

• For d e N*, 7Wd(M) denotes the real d square matrices; 5d(M), 5+(R),5+'*(R), and ^d(R) denote 
respectively the set of symmetric, symmetric positive semidefinite, symmetric positive definite and non 
singular matrices. 

• The set of correlation matrices is denoted by £ci(R): 

e:d(R) = {xe 5+(K),Vl < ^ < d, x^,^ = l} 
We also define £iJ(R) = H Gd(R-), the set of the invertible correlation matrices. 

• For X G AidO^), x'^ , adj(x), det(x), Tr(x) and Rk(a;) are respectively the transpose, the adjugate, the 
determinant, the trace and the rank of x. 

• For X G 5j"(]R), y/x denotes the unique symmetric positive semidefinite matrix such that {y/x)^ — x 

• The identity matrix is denoted by Id- We set for 1 < i,j < d, = {lk=i,i=j)i<kd<d and = e^. 
Last, we define e^*'''^ = + li^jC-'f . 

• For X e we denote by X{i^j^ the value of Xij, so that x = X]i<i<j<d ^{i.j}^^''''"^ ■ We use both 
notations in the paper: notation {xi^j)i<i^j<d is of course more convenient for matrix calculations 
while {xf^i jT^)i<i<j<d is preferred to emphasize that we work on symmetric matrices and that we have 

Xi_j — Xj^i' 

• For Ai, . . . , Ad <E R, diag{Xi, . . . , Ad) G Sd(^) denotes the diagonal matrix such that diag{Xi, . . . , Xd)i i — 

• For X e 5j"(M) such that Xi,i > for all 1 < i < d, we define p(x) G £d(K) by 

{p{x)\,, = l<i,j<d. (5) 

y/Xi^lXjJ 

• For X e 5d(M) and 1 < i < d, we denote by G 5d_i(M) the matrix defined by x^^\ — x^+i^.^ 

and x^ E W^^^ the vector defined by = Xi^k for 1 < fc < i and x], — Xi,k+i for i < fc < d — 1. For 
X £ €d{M.), we have {x - xe^^x)^'^ = a;W - x'{x')^. 
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1 Some properties of MRC processes 

1.1 The infinitesimal generator 

We first calculate the quadratic covariation of MRCd{x, k, c, a). By Lemma [27l we get: 



{d{Xt)^,,,d{Xt)k,i) 



aU^^=k{Xt - XteliXt)jj + l,=i{Xt - Xte'aXt)j,k) 
Wj{lj=k{Xt - Xte'^Xt)^,i + lj=i{Xt - Xte'^Xt)^^k) 



dt 



a,'(l^=fe((Xt)j-i - {Xt)^AXt)^,l) + l^=l({Xt)J,k - (Xt)^^,{Xt)^,k)) 
WAl,^k{{Xt)^.l - {Xt)^,,{Xt)j,i) + lj^l{{Xt\k - 



(6) 



dt. 



We remark in particular that d{{Xt)ij , d{Xt)k,i) = when i,j,k,l are distinct. 

We are now in position to calculate the infinitesimal generator of MRCd{x,K^c,a). The infinitesimal 
generator on A^d(R) is defined by: 

X e e:d(K), L^fix) = lim ^-^^ ^ for / e C'^{Md{R)M) with bounded derivatives. 

t-i.O+ t 

By straightforward calculations, we get from ^ that: 

= X! + '^j)i'^i,3 ~ + \ X! <4{^j,k - Xi^jXi^k)[di,jdi^k + dijdk,i + Oj^idi^k + dj^idk,i\. 



l<i,j<Q 



l<i,j,k<d 



Here, di.j denotes the derivative with respect to the element at the r line and j column. We know 
however that the process that we consider is valued in C Though it is equivalent, it is often 

more convenient to work with the infinitesimal generator on which is defined by: 



X e €diR), Lfix) = lim ^[/(-^f)] /(^) for / e C^Sdi 

t-i.O+ t 



with bounded derivatives, 



Xj^i the value of 



since it eliminates redundant coordinates. For x G 5(j(M), we denote by X[i,jj = Xi.j 
the coordinates and so that x = Y.l<^<j<d^{^,3}(^'d + '^i^j^lH- For / ^ C2(5d(M),M), 

denotes its derivative with respect to x^i^y For x e Md(S.), we set tt{x) = {x + x'^)/2. It is such that 
tt{x) = X for X G 5£i(K), and we have Lf{x) = f o tt{x). By the chain rule, we have for x G 
dijf o tt{x) = {li^j + and we get: 



E - + 1 E "■^(^{jM - X{i,j}xii^k})d{i,j}di,,k} 



i<j<d 



l<j,k<d 
j^i^k^i 



(7) 



Then, we will say that a process {Xt,t > 0) valued in solves the martingale problem of MRCd{x, k, c, a) 

if for any n e N*, < ti < ••• < t„ < s < t, gi,...,g„ G C(5d(M),M), / G C2(5d(M),M) we have: 



E 



f[g,{Xt^)(f{Xt)-f{Xs)-[ Lf{X^)di 



0, and Xq = x 



(8) 



Now, we state simple but interesting properties of mean-reverting correlation processes. Each non-diagonal 
coefficient follows a Wright-Fisher type diffusion and any principal submatrix is also a mean-reverting cor- 
relation process. This result is a direct consequence of the calculus above and the weak uniqueness of the 
SDK ^ obtained in Corollary H 
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Proposition 1 — Let {Xt)t>o ^ MRCd{x,K,c,a). For 1 < i j < d, there is Brownian motion {f3l'-' ,t > 
0) such that 

d(Xt\,, = + nj)(c.^, - {Xt).^,)dt + ^al+a]^\-{Xt)l,dPl''. (9) 

Let I = {ki < ■ ■ ■ < kd'} C {1, . . . ,d} such that 1 < d' < d. For x G Md{R), we define x^ G Md'iJS.) by 
— ^ki,kj for 1 < i,j < d' . We have: 

{Xl)t>o MRCd'ix',^',c',a'). 



1.2 Calculation of moments and the ergodic law 

We first introduce some notations tliat are useful to cliaracterise tlie general form for moments. For every 
X e 5d(M), m e 5d(N), we set: 

X™ = Yl x^^f and \m\ = ^ m{,,,}. 

l<i<i<<i l<i<j<d 

A function / : ^ M is a polynomial function of degree smaller than n G N if there are real numbers 

am such that f{x) = I]|,„|<„ a^a;"', and we define the norm of / by ||/||p = J2\m\<n 

We want to calculate the moments E[X™] of {Xt,t > 0) ~ MRCd{x, k, c, a). Since the diagonal elements 
are equal to 1, we will take m^aj = 0. Let us also remark that for i j such that Ki = Hj ~ 0, we have 
from ([3]) that Oj = 0. Therefore we get {Xt)ij — Xij by ([9]). 



Proposition 2 — Let m G 5(;(N) such that nii^i = for 1 < i < d. Let {Xt)t>o ^ MRCd{x,K,c,a). For 
m G 5d(N), Lx"' = -Kmx"' + f„i{x), with 

d d ^ d d 

i—1 j — 1 2—1 j^k—l 

and 

d d {' '} \ ^ ^ 

4 = 1 j = l 4 = 1 'j,k = l 

is a polynomial function of degree smaller than \m\ — 1. We have 

E [X^] = .T™ exp(-iif„) + exp(-ii^„) / exp(sX„OE[/„(X,)]ds. (10) 

"'0 



Proof : The calculation of Lx™ is straightforward from By using Ito's formula, we get easily that 



dt 



= +E[/„(Xt)], which gives (Hn 



Equation (jlOp allows us to calculate explicitly any moment by induction on \m\. Here are the formula 
for moments of order 1 and 2: 

VI < » 7^ j < d, E [(XO.j] - a;,,,e-*(«'+'^^' + c,,,(l - e-*('='+'=^)), 

and for given 1 < i ^ j < d and 1 < k ^ I < d such that + Kj > and Kk + ki > 0, 

E [(X4)ij(Xt)fe,;] = XijXk,ie"*^'-^-''-' + (Kj + Kj)CijJkd{t) + (Kk + Kl)Ck,ni,jit) 

+a'^ (li=fe7j\i (0 + l4=/7i,fe(i)) + (lj=fe74,i W + li=a»,fc(i)) , 
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where -ftTij.fc,; ^ Hi + kj + Hk + m + af (li^k + li=i) + a] (lj=fc + and 

\ _ *-fi,3,fc,! g — *(Km + ft„) _ gt-R'i,j,fc,! 

Let / be a polynomial function of degree smaller than rt G N. From Proposition [21 L is a linear mapping 
on the polynomial functions of degree smaller than n, and there is a constant C„ > such that ||i/||p < 
Cn||/I|p- the other hand, we have by Ito's formula E[/(Xt)] = f{x) + E[L f {X s)]ds , and by iterating 

E[/(Xt)] = Eto l<L'fix) + /o ^^E[L'=+i/(^s)]rfs. Since ||L7||p < C;,||/|1p, the series converges and we 
have 

OO j 

E[/(xo] = ^^i7(x) (11) 

i=0 

for any polynomial function /. We also remark that the same iterated Ito's formula gives 

V/eC°°(5d(K),M),VfceN*,3C7>0,Vxe CdW, |E[/(XO] - V ^LVWI < (12) 

since L'^'^^f is a bounded function on £d(R)- 

Let us discuss some interesting consequences of Proposition [51 Obviously, we can calculate explicitly in 
the same manner E[X™^ . . . X™"] for < ii < • • • < t„ and toi, . . . , m„ e Sd{N). Therefore, the law of 
(Xjj , . . . , ^t„) is entirely determined and we get the weak uniqueness for the SDE ([2]). 



Corollary 3 — Every solution {Xt,t > 0) to the martingale problem ^ have the same law. 

Proposition [21 allows us to compute the limit limt_j.+oo E[X™] that we note E[X™] by a slight abuse of 
notation. Let us observe that Km > if and only if there is i,j such that Hi + Hj > and m,;.j > 0. We have 

E[X™] = X™ if m e Sd{N) is such that m^.jy > k, = kj = 0, (13) 

E[X™] - E[/„,(X^)]/i^,„ otherwise. 

Thus, Xf converges in law when t -> +oo, and the moments E[X™] are uniquely determined by (|13p with 
an induction on |m|. In addition, if Hi + nj > for any 1 < i, j < d (which means that at most only one 
coefRcient of k is equal to 0), the law of X^o does not depend on the initial condition and is the unique 
invariant law. In this case the ergodic moments of order 1 and 2 are given by: 

E[(Xoo).j] = c,,„ 

wuv \ ^v \ 1 {i^i + 1^] + + Ki)cijCk,i + a'f{li^kCj,i + h^iCj^k) + a'jilj=kCi,i + lj=iCi^k) 

'^[(Xoo)i,j(Xoo)k,l\ — f-; ■ 

1.3 The connection with Wishart processes 

Wishart processes are afRne processes on positive semidefinitc matrices. They have been introduced by 
Bru [5 and solves the following SDE: 

+ {{a + l)a^a + bY^ + Y^b^) ds + (^VY^dW^a + a^dW^^^YF^ , (14) 

where a,b € AidO^) and y € 5j~(]R). Strong uniqueness holds when a > d and y G iSj^'*(R). Weak existence 
and uniqueness holds when a > d — 2. This is in fact very similar to the results that we obtain for mean- 
reverting correlation processes. The parameter a + 1 is called the number of degrees of freedom, and we 
denote by WISdiy, a + l,b, a) the law of (F/, t > 0). 
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Once we have a positive semidefinite matrix y S 5j"(R) such that yi^i > for 1 < i < d, a trivial way to 
construct a correlation matrix is to consider p(y), where p is defined by ([5]). Thus, it is somehow natural 
then to look at the dynamics of p(F/), provided that the diagonal elements of the Wishart process do not 
vanish. In general, this does not lead to an autonomous SDE. However, the particular case where the Wishart 
parameters are a = e\ and 6 = is interesting since it leads to the SDE satisfied by the mean-reverting 
correlation processes, up to a change of time. Obviously, we have a similar property for a = and 6 = by 
a permutation of the iih and the first coordinates. 



Proposition 4 — Let a > max(l,d — 2) and y G ^^"(IR) such that yi^i > for 1 < i < d. Let (F/)t>o ~ 
WISd{y,a + 1,0, e;^). Then, {Yf')i^i = yi,i for 2 < i < d and follows a squared Bessel process of 

dimension a + I and a.s. never vanishes. We set 



Xt = p(r/), m = f jri—ds. 

Jo [J^s )1A 

The function (p is a.s. one-to-one on R-|_ and defines a time-change such that: 

> 0) MRCd{p{y), ^e\,Id,e\). 

In particular, there is a weak solution to MRCd{'p{y),%e-\,Id,e]i). Besides, the processes {X^-i(^t^,t > 0) 
and > 0) are independent. 

Proof : From ([M]), a = e;^ and 6 = 0, we get d(y/),,j = for 2 < i, j < d and 

d d 

= {a + l)dt + 2Y,i^p?)l.k{dWt)k,u = (15) 

k=l k=l 

In particular, ~ 4(y/)i^idi and is a squared Bessel process of dimension a + 1. Since 

a + 1 > 2 it almost surely never vanishes. Thus, {Xt, t > 0) is well defined, and we get: 

diX.h. = -f(^*)M7^+Ef^^^-(^0M^^^) W)m (16) 

By Lemma [5(71 (j){t) is a.s. one-to-one on and we consider the Brownian motion {Wt,t > 0) defined by 
{Wcj,{t))i,j — Jo '^^yf^'^^' have by straightforward calculus 

k=i \ ^y^^' y'(i7-i(t))i,iy 

d{{X^-i(^t))i,i,{^<p-Ht))i,j) = [{X<t,-^(t))i,j - {X<t,-Ht))iA^4>-^(t))i,j]dt^ 

which shows by uniqueness of the solution of the martingale problem (Corollary [3]) that {X^-i(^f),t > 0) 
Mi?C<i(p(2/),fei,/d,ei).^ 

Let us now show the independence. We can check easily that 

= YTpi^iiXtkj - {Xth4Xt)i.,] and d{{Xt),.„ (Yy),,,) = 0. (18) 
['^t 71,1 

We define ^'(y) e 5d(K) for y e Sj{R) such that y^^, > by «'(y)i,j = *(2/)i,i = yia/y/yiJvU and 
^{y)i^j = yij otherwise. By and ([H]), (*(yf),t > 0) solves an SDE on 5<j(M). This SDE has a unique 



7 



weak solution. Indeed, we can check that for any solution {Yt,t > 0) starting from ^{y), {^~^{Yt),t > 
0) ~ WISd{y,o^ + 1,0, e^), which gives our claim since 4* is one-to-one and weak uniqueness holds for 
WISd{y,a + 1,0, e;^) (see 0). Let {Bt,t > 0) denote a real Brownian motion independent of {Wt,t > 0). 
We consider a weak solution to the SDE 

diYt)i^i = (a + l)dt + 2^(Ft)i,idBt, d{Yt\j = for 2 < i, j < d, 

^ (l^tjia ^ \ \/{Yt)isyiA (,^tji,i / 

that starts from Yq = ^'(y). It solves the same martingale problem as ^(i^t), and therefore (^"(14),^ > 
0) (^t, t > 0)- We set = J* jyj^ds. As above, « = 2 . . . , d) solves an SDE driven by 

{Wt,t > 0) and is therefore independent of {{Yt)i^i,t > 0), which gives the desired independence. □ 



Remark 5 — There is a connection between squared-Bessel processes and one- dimensional Wright-Fisher 
diffusions that is similar to Proposition^ Let us consider Z\ = Zi -\- Pit -\- (T-\/Z*'c?i?s, i = 1, 2 two squared 
Bessel processes driven by independent Brownian motions. We assume that /3i, /32, c > and < 2(/3i -I-/32) 
so that Yt = Zl -\- Zf is a squared Bessel processes that never reaches 0. By using ltd calculus, there is a real 
Brownian {Bt,t > 0) motion such that Xt = Zl /Yt satisfies 

dXt = -t- - Xt)^ + aVX,(l-XO^, 

Pi + P2 l^t V 

and we have {dXt,dYt) = 0. Thus, we can use the same argument as in the proof above: we set (f)(t) = 
l/{Ys)ds and get that {X^-i(^tj,t > 0) is a one- dimensional Wright-Fisher diffusion that is independent 
of {Yt,t > 0). This property obviously extends the well known identity between Gamma and Beta laws. This 
kind of change of time have also been considered in the literature by JT^ or flSf for similar but different 
multi- dimensional settings. 



1.4 A remarkable splitting of the infinitesimal generator 

In this section, we present a remarkable splitting for the mean-reverting correlation matrices. This result 
will play a key role in the simulation part. In fact, we have already obtained in 2 very similar properties 
for Wishart processes. Of course, these properties are related through Proposition lU which is illustrated in 
the proof below. 

Theorem 6 — Let a > d — 2. Let L be the generator associated to the MRCd{x, ,Id,a) on and 
Li be the generator associated to MRCd{x, ^e^, Ld, e^), for i G {1, . . . , d}. Then, we have 

d 

L = '^^afLi and yi,j G {1, . . . ,d}, LiLj — LjLi. (19) 

i=l 

Proof : The formula L ~ '^i -^i obvious from The commutativity property can be obtained 

directly by a tedious but simple calculus, which is made in Appendix O Here, we give another proof that 
uses the link between Wishart and Mean- Reverting Correlation processes given by Proposition |4l 

Let Lf denotes the generator of WISdix,a + l, 0, e'J. From pi, we have Lf'Lf = Lf Lf for 1 < i, j < 

d. Let us consider a > max(5,d-2) and x G ^^(K)- We set for i = 1,2 (y/'^,i > 0) M^/S'rf(x, a -f 1, 0, ejj). 
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and we assume that the Brownian motions of their associated SDEs are independent. Since LY LY = L^^ii 
we know from [2 that 1^/'^' y'^"'^* and thus 

E[/(p(y/'^''"))]=E[/(p(r/'^''^))], 

for any polynomial function /. By Proposition |31 p(Y'/'^' ) x]^y-i(^nt))^ where {x]^'^y\i^_^yU > 0) is a 

mean-reverting correlation process independent of (t)^{t) = /„ ds. Since {Y^ i = 1, {Ys ' ' )i i 

follows a squared Bessel of dimension a + 1 starting from 1. Using the independence, we get by (|12l) 

By LemmaEH wehaveE[0i(i)] ^ t + + 0(t^), E[(j)\t)^] = t + 0{t^),E[(l)^{t)^] = 0{t^). Thus, we get: 

E[/(p(r/'^' = /(p(r,2'^)) + iLi/(p(y/'^)) + |[L?/(p(y/'-)) + (3 - a)Li/(p(y,^'-))] + ©(t^). 

Once again, we use Proposition|4]and (IT2|) to get similarly that E[f{■p{Y^'^))] = f (x) +tL2f{x) + ^[L2f (x) + 
(3 — a)L2f{x)] + 0{t^) for any polynomial function /. We finally get: 

E[/(p(r/'^''"))] = /(x) +i2)/(.T) + |[L?/(x) + 2L2Li/(a;) + + (3-a)(Li +L2)/(.t)] +0(^3). 

Similarly, we also have 

E[fipiY,'-''^''^))]=fix)+tiL^+L2)f{x)+*^[Llf{x) + 2L,L2fix) + L^^^^ 

(20) 

and since both expectations are equal, we get LiL2f{x) — L2Lif{x) for any a > max(5,(i — 2). However, 
we can write Li = \ {aLf + L^'), with 

Lf = ^ a;{j,j}9{ij} and Lf^ = ^ {x{]M ^ ^{^■j}^{^M)^{i,j}^{iM- 

l<j<d l<j.k<d 

Thus, we have a^if + aiL^hf + if if ) + if if - a^Lf Lf + aiL^hf + Lf Lf ) + Lf if for any 
a > max(5, d - 2). This gives Lf = Lf , Lf + = Lf + Lf Lf , Lf = , and 

therefore iiL2 = L2L1 holds without restriction on a. □ 

Remark 7 — Let x e €d{M.), {Yt'''',t > 0) - WISd{x,a + 1,0, e^) and L^ its infinitesimal generator. 
Equation ^ and the formula E[f{-p{Y^ = f{x) + tLif{x) + ij[Llf{x) + (3 - a)Lif{x)] + 0{t^) used 
in the proof above lead formally to the following identities for x G £(i(K) and f G C°°{Sd{^), K), 

if (/ o p){x) = LJ{x), [Lffif o p)(x) = Llf{x) + (3 - a)Li/(x), Lf L^ [f o v){x) = LMix), 

that can be checked by basic calculations. 

The property given by Theorem [S] will help us to prove the weak existence of mean-reverting correlation 
processes. It plays also a key role to construct discretization scheme for these diffusions. In fact, it gives a 
simple way to sample the law MRCd{x, ^a^,Id,a;t). Let x G £(i(R)- We construct iteratively: 

. Xl'"" ^ MRCdix, f afei , Id, aie^; i). 
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For 2 < i < d, conditionally to ' , Xl'"' ' ^ MRCd{Xl * ,^a1e\,Id,aie\;t) is sam- 



pled independently according to the distribution of a mean-reverting correlation process at time t with 
parameters {^afe'^j^, Id,aie''^) starting from Xl~^' 



Propositions — Let X^'" " he defined as above. Then, xf'" ' ^ M RCd{x, ^a'^ , Id,a;t). 

Let us notice that MRCd{x, f afe^, Id, aiSd'^ t) MRCd{x, ^e^, Id, e^; aft) and that MRCd{x, f e^, Id, e^; t) 
and MRCd{x, f-e^, /d, e^; i) are the same law up to the permutation of the first and the i-th coordinate. Thus, 



it is sufficient to be able to sample this latter law in order to sample MRCd{x, §a^, Id, a; t) by Proposition[5] 

a„2 



Proof : Let / be a polynomial function and - MRCd{x,^a^ ,Id,a;t). By E[/(Xf 



Ej°lo 7^ Using once again ((II|,E[/(X,'^- ' )] = E[E[/(xf- ' )\X, 
= J2fLo jT^[Ldf{X^ ~ '"' )], and we finally obtain by iterating 



fjl-l hjd jL^ p 



nf{xt'-' )]= -^L{^...L^d'f{x)=Y,-{L, + --- + Ldyf{x)^nf{X^)l 

n Jl - • • ■ Jd- ■ n J • 

since the operators commute. □ 

We can also extend Proposition [5] to the limit laws. More precisely, let us denote by MRCd{x, k, c, a; oo) 
the law characterized by (IT^ . We define similarly for x £ X^^ ^ MRCd{x, f af ej, Id, oiej; oo) and, 

conditionally to X'^'^'-^'^^ , X^ '''"" ~ Mi?Cd(X^^— ""^"^ , f a^ej,, Id, aic'^; oo) for 2 < i < d. We have: 



X^-""^ ^ MRCd{x,^a',Id,a;^). (21) 

To check this we consider {Xt,t > 0) ^ MRCd{x, ^a^ , Id,a) and m S 5d(N) such that m^.i = 0. By 
Proposition [21 E[X™] is a polynomial function of a; that we write E[X/"] — J2m'eSd{N) \m'\<\m\ lm,m'{t)x™ . 
From the convergence in law (1131) . we get that the coefficients 7m, m' (t) go to a limit jm,m' (oo) when t — > +oo, 
and E[X™] = J2\m'\<\rn\lm,m'ioo)x"'' . Similarly, the moment m of MRCd{x,^a1e''^,Id,aie'^;t) can be 
written as Yl,\m'\<\m\ 7m,m'(*)2^"'- We get from Proposition [U 



[Xt \— ^ ^ 7m,md (^)7jTid,md_i (^) ■ • ■ 7m2,mi (0 



mi 



|mi|<---<|mjj|<|m| 

which gives (PT|) by letting t +oo. 

1.5 A link with the multi-allele Wright-Fisher model 

Theorem[S]and Proposition|S]have shown that any law MRCd{x, ^a^,Id, a; t) can be obtained by compo- 
sition with the elementary law MRCd{x, ^,Id,e\;t). By the next proposition, we can go further and focus 
on the case where {xij)2<i.j<d = Id-i- 

Proposition 9 — Let x G £d(K)- Let u G A^d-i(IR) o.nd x G £(i(R) such that x — i \ '^j^fn 'r 



^ u / V 

and {x)2<ij<d = Id-i fLemma \26] gives a construction of such matrices). Then, for a>2, 
MRCd{x,'^e\,Id,e\) J M MRCd{x,^e\,Id,e]^ ( \ J. 
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Proof : Let {Xt,t > 0) ~ Mi?Cd(x, f e^, 7^, ei). We set = J || ^ J J, ^ . Clearly, 

{{Xt)i,j)2<i,j<d = Id-1 and the matrix {{Xt)t,j)2<i,j<d is constant and equal to uu'^ = {xij)2<ij<d- We have 

for2 < i < d, = Efc=2"^-i,fc-i(^0i,fe- By ©, we get (d(Xt)i,fe, ^ [lk=i- iXt)i,kiXt)i,i]dt- 

Therefore, the quadratic variations 

(d d \ 

^ - ^ Ui^i.k-i{Xt)i,kUi-i.i-i{Xt)i,i dt 

k=2 k,l=2 J 

=. {{Xt)^^j ~ {Xt)i4Xt)i,,) dt, 

are by ([6]) the one of MRCd{x, §ej, 7^, e^). This gives the claim by using the weak uniqueness (Corollary [3]). 
□ 

For X e iSd(K) such that {xij)2<i,j<d — Id-i and — 1, we have det(a;) — 1~ ^2^=2 -^i i ^-'^'^ therefore 



a; e 



E^l^^l- (22) 



The process {Xt)t>o ^ MRCd{x, ^,Id,e^;t) is such that {{Xt)ij)2<i,j<d = Id-i- In this case, the only non 
constant elements are on the first row (or column). More precisely, ((^t)i,j)i=2,...,d is a vector process on 
the unit ball in dimension d — 1 such that 

d{{Xt)l,^,{Xt)l,,) = (1,=, - {Xt)l,^{Xt)l,j)dt. 

For i = 1, . . . ,d - 1, we set C*' = (Xf)?,,+i. We have {dC,l,dC,l) = 4(^1^=^ - C,l)dt and the drift of Q is 
(1 — (1 + 2a)Q)dt . Thus, {Ct)i<i<d-i satisfies X^t^i Q — ^ and has the following infinitesimal generator 

^[1 - (1 + 2a)z^]d,^ +2 ^»(1»=J - ^j)9zA, 

i=l l<i.j<d-l 

This is a particular case of the multi-allele Wright-Fisher diffusion (see for example Etheridge [H]), where 
{(,},... Xt ^^ A ~ J2i=i Q) describes population ratios along the time. Similar diffusions have also been 
considered by Gourieroux and Jasiak 14 in a different context. Roughly speaking, {{Xt)i42<i<d can be 
seen as a square-root of a multi-allele Wright-Fisher diffusion that is such that its drift coefficient remains 
linear. 

Also, the identity in law given by FropositionlH] allows us to compute more explicitly the ergodic limit law. 
Let X e e:d(R) such that ixi^j)2<t,j<d ^ Id-i, {Xt)t>a MRCdix, f e^. Id, e]) and (Yf)t>o ~ WISdix, a + 
1,0, e^). We know by [2] that {(Yf)ij)i<i^j<d = Id-i and 



((l^")i,^)i<^<rf {Zt'-' + Y.ixi,, + VtN.f, xi.2 + ViN2, . . .,xi^d + VtNd), 

1=2 



where Ni ~ A/'(0, 1) are independent standard Gaussian variables and Z^^'^ — xi^i + {a + 2 — d)t + 
2 \/ zZ^'^dPu is a Bessel process independent of the Gaussian variables starting from xi^i. By a time 
scaling, we have Zi^'^ tZi^'^^^ , and thus: 

. ^ lav. ( ^ + , • ■ • , ^ + iV,) (iVs, . . . , TVrf) 

(P(j^t )l,i)2<i<d = 



z-^^/* + j:t2C-^ + N.y v^zf + Et2 Nf 
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On the other hand, we know that converges in law when t — )• +cx), and Proposition 2] immediately gives, 
with the help of Lemma l30l that ((X^)i,i)2<i<d /l^" "'^''^ 2 ' simple calculations, we get that 

{{X^)i^i)2<i<d has the following density: 

In particular, we can check that i)2<i<d follows a Dirichlct law, which is known as the ergodic limit 

of multi-allele Wright-Fisher models. Last, let us mention that we can get an explicit but cumbersome 
expression of the density of the law MRCd{x, ^a^, /^j, a; c«) by combining (I?!]) . Proposition [S] and 

2 Existence and uniqueness results for MRC processes 

In this section we show weak and strong existence results for the SDE ©, respectively under assump- 
tions (O and (U]). These assumptions are of the same nature as the one known for Wishart processes. To 
prove the strong existence and uniqueness, we make assumptions on the coefficients that ensures that Xt 
remains in the set of the invertible correlation matrices where the coefficients are locally Lipschitz. This is 
similar to the proof given by Bru 5 for Wishart processes. Then, we prove the weak existence by introducing 
a sequence of processes defined on £rf(]R), which is tight such that any subsequence limit solves the martingale 
problem (|H]). Next, we extend our existence results when the parameters are no longer constant. Last, we 
exhibit some change of probability that preserves the global dynamics of our Mean-Reverting Correlation 
processes. 

2.1 Strong existence and uniqueness 

Theorem 10 — Let x £ £^(M). We assume that ^ holds. Then, there is a unique strong solution of the 
SDE © that is such that yt > 0,Xt G €*j^{R). 

Proof : By Lemma [H we have (^x - xe^-x)'"' v^xW - and xl"! - e (R) when 

X e e:;j(R). For x e such that xl"! - a;"(a;")^ G Sjl\{M), we define /"(a;) € Sj{B) by (/"(a;))„j = 

for 1 < j < d and (/"(x))'"' = -y/ — a;"(a;")^. The function /" is well defined on an open set of iSd(M) 
that includes and is such that /"(x) = ^JlE^^lL^^ for x e Since the square-root of a positive 

semi-definite matrix is locally Lipschitz on the positive definite matrix set, we get that the SDE 

/.f d „t 

Xt = x+ {k{c - Xs) + {c~ Xs)n) ds + Y,an (/"(X,)dW^,eS + eSdV7j/"(X,)) , 
■^0 n=i -^0 

has a unique strong solution for < i < r, where 

r = inf{i >0,Xt<^ 5+'*(R) or 3i G {1, . . . ,4,^4^ - X^Xlf ^ S+l\{R)}, inf = +oo. 
For 1 < i < d, we have (/" (X,)dVK,e3),M = /"(^s)„j(dVK,)j,„ = and then: 

d{Xt)^^^ = 2k,41 - {Xt)^,i)dt, 

which immediately gives (Xt)»,j = 1 for < i < r. Thus, Xt G Q(R) for < t < r and t = inf{t > 0, ^ 
£^(R)} by Lemma and the process Xt is solution of ^ up to time r. We set Yt = log(det(Xt)) +Tr(2K — 
a^)^. By Lemma [28l we have 

Yt = Yo+ f Ti[X-\Kc + cK--da'^)]ds + 2 f JTr[a^X^^ ~ Id)]dl3s 
Jo Jo 

> Yo + 2 f JTr[a^X^' - Id)Ws, 
Jo 
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since nc + cn — da? £ by Assumption Now, we use the McKean argument exactly like Bru [5] did 

for Wishart processes: on {r < oo}, Yt —cx), and the local martingale /* i/Tr[a2(Xj^^ — /d)]d/?s — > — oo, 

t — >T V t — 

which is almost surely not possible. We deduce that r = +00, a.s. □ 



2.2 Weak existence and uniqueness 

The weak uniqueness has already been obtained in Proposition [31 and we provide in this section a 
constructive proof of a weak solution to the SDE In the case d = 2, this result is already well-known. 
In fact, by Proposition [Tl the associated martingale problem is the one of a one-dimensional Wright-Fisher 
process. For this SDE, strong (and therefore weak) existence and uniqueness holds since the diffusion 
coefRcient is l/2-H61derian. 

Thus, we can assume without loss of generality that d > 3. The first step is to focus on the existence 
when a — diag{ai, . . . , Ud) G 5j"(R), a > d — 2, n — and c = Id- By Proposition 21 we know that weak 
existence holds for MRCd{x, ^e^, Id, e]^), and thus for MRCd{x, f ejj. Id, Oie^) for « = 1, . . . , d and > 0, 
by using a permutation of the coordinates and a linear time-scaling. Therefore, by using Proposition |8l the 
distribution MRCd{x, ^a'^,Id, a; t) is also well-defined on £d(R) for any t>Q. Let T > be a time-horizon, 
TV G N*, and tf iT/N. We define {X^ .t £ [Q,T]) as follows. 

• We set = X. 

• For i — Q, . . . ,N ~1, X!^ is sampled according to the law MRCdiXj^ , ^a^. Id, a; TIN), conditionally 
toX^^. 

" roi L t [tj , tj^iJ' ~ T/N t" ~'~ T/N t™|_i ~ t" T/Af ^ tf^i )■ 

The process e [0,T]) is continuous and such that almost surely, Vt G [0,T],Xl^ G £d(R). We endow 

/ ^ \l/4 „ 

the set of matrices with the norm ||x|| = j=i^i,j) ■ The sequence of processes {X^,t G [0,T])jv>i 
satisfies the following Kolmogorov tightness criterion. 

Lemma 11 — Under the assumptions above, there is a constant K > such that: 

yO<s<t<T, m\X^ - X^W^] < K{t - sf. (24) 



Proof : We first consider the case s — t^ and t — tf' for some < k < I < N. Then, by Proposition [51 
we know that conditionally on X!^, X!^ follows the law of MRCdiX!^ , ^a^ , Id, a). In particular, each 

element {X^)i,j follows the marginal law of a one-dimensional Wright-Fisher process with parameters given 
by equation ([9]). Thus, by Proposition [29] there is a constant still denoted by X > such that for any 
l<i,j<d, E[((X^).j - (X^),,j)1 < K{t^ ~ t^f, and therefore 

I k 

n\^-x^^r]<Kd'it? -t^f- 

Let us consider now < s < t < T. If there exists < fc < — 1, such that s,i G [tk ^^k+l\^ then 

E[||Af - Af ||4] = (^^^^E[||X/X - X^mW^] < Kd^{s - tf. Otherwise, there are fc < / such that 

t^,-T/N<s<t«,<tf<t<ti^ +T/N, and E[|| Af - Iff] < Kd^{t^ sf + (t - t^' + (if - 1^?] < 
K'(t - s)^ for some constant K' > Q. □ 
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The sequence {X^,t S [0,T])Ar>i is tight, and we will show that any limit of subsequence solves the 
martingale problem ([8]) . More precisely, we will show that for any nEN*,0<ti<---<tn<s<t<T, 
...,(?„ e C(5d(M),R), / e C°°(5d(R),R) we have: 



lim E 

Af-H-oo 



f[ g^iX^J {j{X^) /(If) - LfiX^)dr 



= 0. 



We set k'^{s) and {t) the indices such that t^j-^-) 



T/N < s <t 



N 



and t 



N 

/"(t) 



< t < t 



Clearly, / is Lipschitz and Lf is bounded on £ci(M). It is therefore sufficient to show that 



N 

/"(t) 



lim I 

N-f+oo 



n.9.(^f) /(^,\ )~/(^: 

We decompose the expectation as the sum of 

n 



N 



i« (t) 



Lf{X^)du 



^ 0. 



(25) 
T/N. 

(26) 



E 



{LfiXl )-LJ{X^:))du 



-E 



■i=l \j=k«{s) 

(27) 

To get that the first expectation goes to 0, we claim that: 



|/?(u,X^^)-/3(tf„(„),X,^ )\du 



(28) 



when /3 : {t,x) e [0,T] x £d(K) — >■ K is continuous. This formulation will be reused later on. By Lemma [TTl 
((28)) holds when j3 is Lipschitz with respect to (t,x). If /? is not Lipschitz, we can still approximate it 
uniformly on the compact set [0,T] x £(i(R) by using for example the Stone- Weierstrass theorem, which 
gives (|28l) . 

On the other hand, we know by ([T2|) that the second expectation goes to 0. To be precise, (fT2l) has been 



obtained by using Ito's formula while we do not know yet at this stage that the process MRCd{x, ^a'^, Id, a) 
exists. It is nevertheless true: p2|) holds for MRCd{x, ^a^e^,/d,e^) since this process is already known to 



be well defined, and we get by using Proposition [5] and Proposition [TH] that 3K > 0, |/(X „ ) — f{X 

{T/N)Lf{X^)\ < K/N"^. Thus, {X^ ,t e [0,T])jv>i converges in law to a solution of the martingale 

problem ([5]). This concludes the existence of MRCd{x, ^a?,Id, a). 

Now, we are in position to show the existence of MRCd{x, K,c,a) under Assumption ([3]). Wc denote by 
^(t, x) the solution to the linear ODE: 



^'(t, x) — k{c ~ x) + {c — x)k 



d-2 



a^ild -~x) + {Id - x)a^], ^iO,x) ^ x e Gidi 



(29) 



By Lemma mi we know that Vt > 0, ^'(t, x) E £d(M). It is also easy to check that: 

3K > o,Vx e e:d(R), U{t,x) - x\\ < Kt. 
Now, we define {X^ ,t £ [0,T]) as follows. 

• We set X^ ^xE €d{^). 

• For i = 0, . . . , iV- 1, Xf is sampled according to MRCd(S,(T/N, Xf ), ^a^, 7^, a; T/N), condition- 
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ally to X^. More precisely, we denote by [X^ ,t £ [tf ,tf_^i]) a solution to 

Xr ^ ^T/N, X^.) + ^ f [a\h - X^) + {h - du 

+ E «n r (^X^-X^e-X^dWue-, + e^W^ ^ X^ - X^ e^X, 

n=l ^ 

and we set Xff, = Xfi, . 

.Forte [t^,tf[,l Xr = X^. + '-^{%^ - 

We proceed similarly and show that the Kolmogorov criterion ([M)) holds for {X^ ,t e [0, r])Ar>i. As already 
shown in Lemma [TTl it is sufficient to check that this criterion holds for s ^ < t ^ tf^ . We have 

WX^.-^r = \\Y.X^^^-^{T/N,X,l)+aT/N,X,l)~X^.f 

< 23 ( wyx^. -x^.\\^ + ii-ky 

\ ^ ^ J + l 3 

\ J=k 

Since [X^ ,t G [0,T]) is valued in the compact set we get easily by using Burkholder-Davis-Gundy 

inequality that E[|| V'-i Xfi - X^^ ||4i < ^(^^ _ ^^ven E[||X^ - Xfi < K{ti - tkf for some 

constant K > Q that does not depend on N . 

Thus, {X^,t e [0,T])jv>i satisfies the Kolmogorov criterion and is tight. It remains to show that any 
subsequence converges in law to the solution of the martingale problem (|S]). We proceed as before and reuse 
the same notations. From (l?7|) . it is sufficient to show that 

3K > 0, \f{X^) - fiX^.) - {T/N)LfiX^.)\ < K/N\ 

Once again, we cannot directly use (|12l) since we do not know at this stage that MRCd{x,K,c,a) exists. 
We have L = + L, where is the operator associated to £,{t,x) and L is the infinitesimal generator of 
MRCdix, ^a^,Id,a). We have: 3K > 0,Va; e CdW, 1/(^(^,2:)) - f{x) - tL^f{x)\ < Kt^ , and 1^ holds 
for L. By Proposition [H we get: 3K > 0,Vx e €d{M.), \f{^{t,x)) - f{x) - tf{x)\ < Kt^, which gives ^ 
and concludes the proof of the weak existence. 

Theorem 12 — Under assumption ([3]), there is a unique weak solution {Xt,t > 0) to SDE ^ such that 

P(Vi>0,Xte£<j(K)) = l. 



Remark 13 — Assumption ^ has only be used in the proof of Theorem ] 12\ to ensure that ^ defined by (|29p 
satisfies 

vt > 0, X e €di^), ^{t, x) e e:rf(K). (30) 

As pointed by Remark \21[ this is a sufficient but not necessary condition. In fact, a weak solution of ([2| 
exists under pO]) . which is more general but less tractable condition than (jSj. 
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2.3 Extension to non-constant coefficients 

In this paragraph, we consider the SDE ^ with time and space dependent coefficients: 

Xt = x+ f [K{s,Xs){c{s,X,)-Xs) + {c{s,Xs)-Xs)K{s,X,)]ds (31) 
Jo 

+ Va„(s,X,) / ^Xs - Xse^XsdW.e"^ + e^^dWjyJX^ - X.e^X 

where K(t, x), c{t, x) and a(t, x) are measurable functions such that for any t > and a; € £^(1^), nit, x) and 
a{t,x) are nonnegative diagonal matrices and c{t,x) € Then, under the following assumption 

VT>0, sup <^, yt e [0,T],3K > 0,\\f{t,x) - f{t,y)\\ < K\\x ~ y\\ for / e {«, c, a}, 

te[o,T] 

yt>0,x& CdW, Kit,x)c{t,x) +c{t,x)K{t,x) - da^{t,x) e 5j"(M) and G C^W, (32) 

strong existence and uniqueness holds for pip . To get this result, we observe that p(a;) is Lipschitz on 
{x e 5+(M) s.t. \/l<i<d, 1/2 < xu < 2}. Therefore, the SDE Xt = x + /J {k{s,p{Xs))[c{s,p{Xs)) - X,,] 
+ [c(s,p(X,)) ~ XM-^,P{Xs)))ds + ELi/o «n(s>P(^s)) (/"(^.)dM^.e2 + eSdW^Jr(^.)) has a unique 
solution up to time r = inf{t >0,Xt^ S^'*{R) or 3i G {1, . . . Xl{XlY ^ 5+4 (M) or (Xf)»,j ^ 

[1/2, 2]}, and we proceed then exactly as for the proof of Theorem [TOl 
Also, weak existence holds for (|3T|) if we assume that: 



K.{t, a;), c(t, x), a(t, a;) are continuous on M_|_ x £d(M) 

Vt > 0, X e CdW, a;)c(t, x) + c(i, a;)K(t, x) - (d - 2)a^(t, x) G 5j'(R). (33) 

To get this result, we proceed as in Section ED and define {X^ ,t & [0,T]) as follows. 

• We set X^ = X. 

• For i = 0, . . . , iV - 1, we denote by {X^ ,t e [if, if_i]) a solution to 



* n{tf,X^.){c{tf,X^^) - X^) + (c(tf ,1,^) - X^Ht^,X^.) 



du 



+ J2<^nitf,X^.) j[ UxN - X^e-^X^dW^e:, + e^^dW^ ^ X^ - ^.^e^JX,^ 
and we set = . 

• For t e [t; Xt = + jy^(X^JV ^ - -'^(n)- 

We can check that {X^ , t e [0, T]) satisfies the Kolmogorov criterion and is tight. To obtain we proceed 
as in Section [2^ More precisely, let us denote for u G [0,T] L„ the infinitesimal generator of (|3T|) . and in 
the infinitesimal generator with frozen coefficient at {ti,Xj^) when u € [tf^ ,tf^_^^). In (P?]). the first term 

E 



t 

N 

before. 



thanks to ([28|) . and the second term goes to as 



To sum up, it is rather easy to extend our results of strong existence and uniqueness, and weak existence 
when the coefficients are not constant. However, we can no longer get explicit formulas for the moments in 
this case. Thus, if the coefficients satisfy p3p but not p2p . the weak uniqueness remains an open question, 
which is beyond the scope of this paper. 
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2.4 A Girsanov Theorem 

In this section, we will use an alternative writing of the SDE ©. In fact, by Lemma (|27p . the SDE 



Jo Ja 



(34) 



is associated to the same martingale problem as MRCdix, K,c^a) for any functions /i„ : 5d(M) — > Aidi^) 
such that hn{x)hn{x)^ = x — xe^x for x € In this paper, we have arbitrarily decided to take the 

symmetric version hn{x) = y/x — xe^x. Obviously, other choices are possible. An interesting choice is the 
following one: 

x G S+[n),K{x) = V^Y^/d- V^e^V^ = V^{Id - V^C^^), (35) 

where the second equality comes from Lemma [24l Obviously, our weak existence and uniqueness results (The- 
orem [121) applies to (p4|) since ([2]) and p4p solve the same martingale problem. However, we have to show 
again that strong uniqueness holds for (IMl) under Assumption (U) and x E £5(R)- The proof is in fact very 
similar to Theorem [TUl We know that there is one strong solution to Xt — x + J* {k{c — Xs) + (c — Xs)k) ds 
+ En=ia™/o (\/^(/d- V^e^V^)dW,e2 + e2(/d-%/^e'iV^)v^(iM^J) up to time r = inf{t > 
0, Xt ^ 5j'(K)}. On t e [0, r), there are real Brownian motions PI such that 



diXt),,^ = 2k,(1 - {Xt),,i)dt + 2a,(l - {Xt)^^,)^/iXt)~,dpi, 

which gives {Xt)i^i = 1 by strong uniqueness of this SDE. We then conclude as in the proof of Theorem [TUl 

and get in particular that Xt E for t>0. 

We consider now a solution to ([34|. and a progressively measurable process {Hs)s>o, valued in Alrf(M), 

such that „ / /■* . T ^ 1 /■* . T ^ ^ 

=exp(^y Tr{HjdW,)~-j Ti{HjH,)ds 

is a martingale. For a given time horizon T > 0, we denote by Q the probability measure, if it exists, defined 

where (J^f )t>o is the natural filtration of the process {Xt)t>o- Then, = Wt — Jq Hgds is a dx d Brownian 
matrix under Q, and the process {Xt)t>o satisfies 



(36) 



Xt = x+ {k{c- Xs) + {c- Xs)K)ds 



(38) 



+ 



' a 



4=1 



XsSdy Xs 



XsGdV Xs 



Hse\ + e'dHi 



ew + e 



Id 



XseWXs 



ds 



'X, 



We present now changes of probability such that (Xt , i > 0) is also a mean- reverting correlation process 
under O. 



Proposition 14 — We assume ([3|). We consider {Xt,t > 0) MRCd{x,K,c,a) and take Ht = \/XtX, 
with A = diag(Xi, . . . , Ad) € iSd(R). Then, p6p is a martingale and {Xt, t > 0) ^ MRCd{x, k, c, a) under Q. 

Proof : Since the process {Xt,t > 0) is bounded, ([5S)) is clearly a martingale. For y e £(i(R): e^ye^ — 
and we have {y/y{Id - ^<^\^)) y^^c^d ^ - y<^dy)'^d = 0' which gives the result by (|3H1)- □ 
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Proposition 15 — Let x E £^(R). We consider {Xt,t > 0) ^ MRCd{x, ,a) and assume that K^,c^,a 
satisfy (|4|). Let (? G £ci(M) and he a real diagonal matrix such that Ui = nf = and k^c^ +c^k^ + 

K^c^ + c^K^ - da^ e 5+(R). We set: 

A = diag{Xi, . . . , A^) with Xi — \ il ^ ^ ^j^fj. — (^^\fx'^^^ (? X 

I otherwise 

This defines with (|36l) and ([37]) a change of probability such that 

{Xt,t>0) ^ MRCdix,K,c,a) under Q, 
where k — diag{Ki, . . . , Kd) G 5^(M) and c G are defined as in Lemma [2R 

Proof : We have ajy/y{Id — y/ye'^y/y)^^^c'^Xe]i — /if (c^e^ — yej^c^e^) = k|(c^ — 2/)e^, which gives the claim 
by psp . provided that E[£^] — 1 for any T > 0. We prove now this martingale property with an argument 
already used in Rydberg [26] and Cheridito, Filipovic, and Yor ([29], Theorem 2.4). 

Let {Xt,t > 0) (resp. {Xt,t > 0)) be a strong solution to ([34]) with parameters k^, c^, a (resp. k, c, a) 
and Brownian motion {Wt,t > 0). For e > 0, we define: 

= inf{t > 0,det(Xt) < e}, HI = lr^>t{^/x'ty^c^X. 

We have limg_^o+ = +oo, a-s. and therefore 

E[f|^] = limE[f|?'l^c>T]. 

On the other hand, we have E[f|?l^e>T] = E[£|f' l^e>T]. We clearly have E[£:|?'] = 1 and = 
Wt — Hgds is a Brownian motion under = £"1?^. Let {Xf,t G [0,T]) be the strong solution to ([34]) 
with the Brownian motion and parameters n, c, a. By construction, X^ = Xt for < t < T A t'^ 
and thus l^e>T = lf->T, where = inf{t > 0,det(Xi^) < e}. We deduce that E[f|f"l^.>T] = Q^f^ > 
T) = P(inf{i > 0,det(Xt) < e} > T) ^ 1, since KC+CK-da^ = Kici+ciKi+K2c2 + c2/«2_da2 ^ 5+ (K). □ 

Let us assume now that > for any 1 < i < d. A consequence of Proposition [TS] is that the probability 
measures induced by MRCd{x,n,c,a) and M RCd{x, k' , c\ a) are equivalent as soon as @ holds for K,c^a 
and K\c\a. By transitivity, it is in fact sufficient to check this for k' = and c' — Ld- By Lemma I22[ 
there is a diagonal nonnegative matrix k and c G such that kc + ck = kc + ck — da^. We get then the 

probability equivalence by using twice Proposition [15] with = |a^, = Ld, = k, = c and = k, 

I 2 ~ 2 ~ 

C C, hyr Kr, C C. 

3 Second order discretization schemes for MRC processes 

In the previous sections, we focused on the existence of Mean-Reverting Correlation processes ([2]) and 
some of their mathematical properties. From a practical perspective, it is also very important to be able 
to sample such processes. By sampling, we mean here that we have an algorithm to generate the process 
on a given time-grid. Through this section, we will consider for sake of simplicity a regular time grid 
tf = iT/N, i = 0, . . . , N ioi a given time horizon T > 0. Despite our investigations, the sampling of the 
exact distribution does not seem trivial, and we will focus on discretization schemes. Anyway, discretization 
schemes are in practice equally or more efficient than exact sampling, at least in the case of square-root 
diffusions such as Cox-Ingersoll-Ross process and Wishart process (see respectively |4] and [2]). First, let us 
say that usual schemes such as Euler-Maruyama fail to be defined for ([2]) as well as for other square-root 
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diffusions. Indeed, tliis scfieme is given by 



X. 



N 



X^l + {k{c-X^I) + {c~X._ 



N ' 



' ' N 



(39) 



+ 



vN 



X 



vN 



X 



N . 



Thus, even if Xfi e Xfi^ can no longer be in and the matrix square- root can no longer be 

defined at the next time-step. A possible correction is to consider the following modification of the Euler 
scheme: 

X,^^^=p((X,«J+), (40) 

where X^^ denotes the right hand side of Here, G 5^(R) is defined for x £ iSd(M) as the unique 

symmetric semidefinitc matrix that shares the same eigenvectors as x, but the eigenvalues are the positive 
part of the one of x. Namely, x^ ~ odiag{X^ , . . . , Aj|)o for x G 5(i(M) such that x — odiag{Xi, . . . , Xd)o 
where o is an orthogonal matrix. Let us check that this scheme is well defined if we start from X!^ e £<i(R)- 

By Lemma [23l the square-roots are well defined, we have (Xjjv)i.i ~ 1 and thus (Xjjv)^^ > 1 and p((Xjjv)+) 
is well defined. By induction, this modified Euler scheme is always defined and takes values in the set 
of correlation matrices. However, as we will see in the numerical experiments, it is time-consuming and 
converges rather slowly. 

In this section, we present discretization schemes that are obtained by composition, thanks to a splitting 
of the infinitesimal generator. This technique has already been used for square-root type diffusions such as 
the Cox- IngersoU- Ross model [4] and Wishart processes [2], leading to accurate schemes. The strength of 
this approach is that we can, by an ad-hoc splitting of the operator, decompose the sampling of the whole 
diffusion into pieces that are more tractable and that we can simulate by preserving the domain (here, the 
set of correlation matrices). Besides, it is really easy to analyze the weak error of these schemes. 



3.1 Some results on the weak error of discretization schemes 

We present now the main results on the splitting technique that can be found in [3] and [5] for the 
framework of Affine diffusions. Here, we have in addition further simplifications that comes from the fact 
that the domain that we consider E) C is compact (typically £ci(R) or D = {a; £ ^"^^^ ^1} in 
Appendix p. For 7 e N'^, we set d^f = d}^^ ...d^^ and I7I = ELi7»- We denote by C°°(D) the set of 
infinitely differentiable functions on D and say that that (C-y)^gNC is a good sequence for / S C°°(D) if we 
have maxxeo\d'yf{x)\ < C^. A differential operator Lf{x) = X]o<|7|<2 '^7(^)^7-/'(^) satisfies the required 
assumption if we have a-y G C°°(D) for any 7. This property if of course satisfied by the infinitesimal 
generator ([7]) of MRCd{x,K,c,a) since the functions are either affine or polynomial functions of second 
degree. Since we are considering Markovian processes on B, we will by a slight abuse of notation represent a 
discretization scheme by a probability measure Px{t){dz) on D that describes the law of the scheme starting 
from a; G B with a time step t > Q. Also, we denote by Xf a random variable that follows this law. Then, 
the discretization scheme on the full time grid (if , i = 0, . . . , N) will be obtained by: 

• xfi=xeO, 

• conditionally to X!^, Xj^ is sampled according to the probability law p<>n {T/N)(dz), and we write 



with a slight abuse of notation X 




A discretization scheme is said to be a potential v-th order scheme for the operator L if for a sequence 
(C-y)^gpj< G (R_i_)^'' , there are constants C, 77 > such that for any function / G C°°(D) that admits (C-y)^gfjc 
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as a good sequence, we have: 



yt e [0,7]), X e 



k=l 



< Ct 



u+1 



(41) 



This is the main assumption that a discretization scheme should satisfy to get a weak error of order v. This 
is precised by the foUowing theorem given in [4] that rehes on the idea developed by Talay and Tubaro [28] 
for the Euler-Maruyama scheme. 



Theorem 16 — Let L he an operator satisfying the required assumptions on a compact domain D. We 
assume that: 

L is a potential weak vth-order scheme for L, 

2. / : D — > M is a function such that u{t,x) = E[/(J("^_j)] is defined and C°° on [0,T] x and solves 
yt e [0, T] , Vx e B, dtu{t, x) = -Lu{t, x) . 

Then, there is K > 0, No e N, such that |E[/(X^)] - E\f(XS,)]\ < K/N" for N > Nq. 

The mathematical analysis of the Cauchy problem for Mean-Reverting Correlation processes is beyond 
the scope of this paper. This issue has recently been addressed for the case of one-dimcnsional Wright-Fisher 
processes by Epstein and Mazzeo [T0|, and Chen and Stroock [7] for the absorbing boundary case. In this 
setting, Epstein and Mazzeo have shown that u{t,x) is smooth for / S C°°([0,1]). However, since we have 
an explicit formula for the moments (jlO[) . we obtain easily that for any polynomial function /, the second 
point of Theorem 1161 is satisfied. By the Stone- Weierstrass theorem, we can approximate for the supremum 
norm any continuous function by a polynomial function and get the following interesting corollary. 

Corollary 17 — Let X^ be potential weak vth-order scheme for MRCd{x, K,c,a). Let f be a continuous 
function on Then, 

Ve > 0,3i^ > 0, \E[f{X^^)]-E[f{X^)]\ <e + K/N''. 

Let us now focus on the first assumption of Theorem 1161 The property of being a potential weak order 
scheme is easy to handle by using scheme composition. This technique is well known in the literature and 
dates back to Strang [27] the field of ODEs. In our framework, we recall results that are stated in [4]. 

Proposition 18 — Let Li,L2 be the generators of SDEs defined on D that satisfies the required assumption 
on D. Let Xl'^ and Xf'^ denote respectively two potential weak vth-order schemes on D for Li and L2. 

1. The scheme xf''^* is a potential weak first order discretization scheme for Li + L2- Besides, if 
L1L2 = L2L1, this is a potential weak vth-order scheme for Li + L2- 

2. Let B be an independent Bernoulli variable of parameter 1/2. If v>2, 

(a) BX;'"^' + (1 - and (b) X^^^* 

are potential weak second order schemes for Li -t- L2- 

Here, the composition X^^ ^ means that we first use the scheme 1 with time step ti and then, conditionally 
to Xl^^ , we sample the scheme 2 with initial value X^^^ and time step ^2- 
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3.2 A second-order scheme for MRC processes 

First, we split the infinitesimal generator of MRC'd{x, n, c, a) as the sum 

L = L« + Z, 

where L is the infinitesimal generator of MRCd{x, ^^a^,Id,a) and is the operator associated to £,{t,x) 
given by Obviously, the ODE (1^9)) can be solved explicitly and we have to focus on the sampling of 

MRCd{x, ^-^a^jldya). We use now Theorem |6] and consider the splitting 

d 

L = ^afL,, 

1=1 

where Li is the infinitesimal generator of MRCd{x, ^^e^, /d, e^). We claim now that it is sufficient to have 
a potential second order scheme for MRCd{x, ^^e^. Id, e;^) in order to get a potential second order scheme 
for MRCd{x, K, c, a). Indeed, if we have such a scheme, we also get by a permutation of the coordinates a 
potential second order scheme X^''^ for MRCd{x, ^^e^, Id, ejj). Then, by time-scaling, X^'j^^ is a potential 

second order scheme for MRCd{x, ^^af e^. Id, CLiC^^). Thanks to the commutativity, we get by Proposition[T51 
that X 2 ■■ is a potential second order scheme for L. Last, still by using Proposition [TS] we obtain that 

S,{t/2, X^2f ) is a potential second order scheme for MRCd{x, n, c, a). (42) 

Now, we focus on getting a second order scheme for MRCd{x, ^^e^, Id,e^). It is possible to construct 
such a scheme by using an ad-hoc splitting of the infinitesimal generator. This is made in Appendix |DJ Here, 
we achieve this task by using the connection between Wishart and MRC process and the existing scheme 
for Wishart processes. In Ahdida and Alfonsi j2], we have obtained a potential second order scheme Y^'^ 
for WISd{x,d — 1,0,6;!;). Besides, this scheme is constructed with discrete random variables, and we can 
check that there is a constant A' > such that for any 1 < i < d, \{Yf. — 1| < AT-s/i holds almost surely 
(we even have {Y^^•'')^,^ = 1 for 2 < i < d). Therefore, we have 1/2 < (F/''')^,^ < 3/2 for t < l/iAK"^). Let 
/ e C°°(e:d(M)). Then f{p{y)) is C°° with bounded derivatives on {y e 5+(R),l/2 < y,^i < 3/2}. Since 
y/'^ is a potential second order scheme, it comes that there are constants C, 77 > that only depend on a 
good sequence of / such that 



Vie (0,77), 



< Ct^, (43) 



where LY is the generator of WISdix,d — 1,0, e;^). Thanks to Remark [3 we get that there are constants 
C, rj depending only on a good sequence of / such that 



Vte (0,7/), 



2 



E[/(p(y/'-))] - !{x) - ( t + (5 - d)M hj{x) - ^-{hff{x) 



< Ct^. (44) 



In particular, p(Y^ '^) is a potential first order scheme for Li and even a second order scheme when d = 5. 
We can improve this by taking a simple time-change. We set: 

[t-(5 -d)^if d>5 

<P[l^) - S -l + ^l+2(5-d)t 

[ ^ — otherwise, 

so that in both cases, 0(t) — t— {b — d)^ + 0{t'^). Then, we have that there are constants C, rj still depending 



only on a good sequence of / such that Vi G (0,77), ¥,[/ {p{Y^^^^))] — f{x) — tLif{x) — ^{LiY f{x) 
and therefore 



< Ci3, 



P(^im) ^ potential second order scheme for MRCd{x, — - — e^. Id, e;^). (45) 
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3.3 A faster second-order scheme for MRC processes under Assumption ( H6i) 

We would like to discuss on the time complexity of the scheme given by ([^^ and (|^5|) with respect to the 
dimension d. The second order scheme given in Ahdida and Alfonsi [5] for WISd{x, d— 1, 0, e^) requires 0{d^) 
operations. Since it is used d times in (I12t to generate a sample, the overall complexity is in O(d^). In the 
same manner, the second order given in Appendix [D] requires 0{d'^) operations. However, it is possible to 
get a faster second order scheme with complexity O(d^) if we make the following assumption: 

ai = ■ ■ ■ ~ ad (i.e. a — aild) and kc + ck— {d — l)a^ G iSj"(R). (46) 

This latter assumption is stronger than ([3]) but weaker than (|4]), which respectively ensures weak and strong 
solutions to the SDE. Under (jm, we can check by Lemma that 

C{t, x) = n{c - x) + (c - x)k - ^-^[a^{Id -x) + (Id - x)a^], ({0, x) = x e (47) 

takes values in £rf(R). Then, we split the infinitesimal generator of MRCd{x, k, c, a) as the sum 

L = + a\L, 

where is the operator associated to the ODE and L is the infinitesimal generator of MRCd{x, ^^Id, Id,Id)- 
In [2], it is given a second order scheme 1^^ for WISd{x,d,0, Id) that has a time-complexity in 0{d'^). We 
then consider / S C°°{<td{^)) and get by using the same arguments as before that there are constants 
C, 77 > depending only on a good sequence of / such that 



Vie (0,7?), 



2 



E[f{p{Yn)] - fix) - tL'^if o p)(x) - ^-{L^ j'if o p)(x) 



where is the infinitesimal generator of WISd{x,d,0, Id)- Thanks to Remark[71 we get that 



ViG (0,7?), 



E[/(p(y,-))] - fix) -lt + i4-d)- Lf{x) - -L'fix) 



In particular, piYf) is a first order scheme for MRCdix, ^^Id, Id, Id) and by Proposition [151 

(it, piY^2f)) is a potential first order scheme for MRCdix, k, c, a). (48) 

As before, we can improve this by using the following time-change: ipit) = i — (4 — d)^ if d > 4 and 
^^^^ ^ -i+^i+2(4^d )t_ Q^j^gj.^igg^ gQ ^^^^ i4:~d)^+ Oit^) in both cases. We get that piY^^^t)) is 

a potential second order scheme for MRCdix, ^^Id,Id,Id)- Then, we obtain that 

C(i/2, P(^a'2^^(f)^'')) is a potential second order scheme for MRCdix, k, c, a) (49) 

by using Proposition [THl Its time complexity is in Oid'^). 



3.4 Numerical experiments on the discretization schemes 

In this part, we discuss briefiy the time needed by the different schemes presented in the paper. We 
also illustrate the weak convergence of the schemes to check that it is in accordance with Corollary [iTl In 
Table [1] we have indicated the time required to sample 10® scenarios for different time-grids in dimension 
d = 3 and d = 10. These times have been obtained with a 2.50 GHz CPU computer. As expected, the 
modified Euler scheme given by (1401) is the most time consuming. This is mainly due to the computation 
of the matrix square-roots that require several diagonalizations. Between the second order schemes that are 
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d = 3 


d= 10 




19 


224 


2nd Qj.jgj. 


65 


1677 


2nd Qj.^gj. "direct" 


90 


3105 


1** order "fast" 


19 


224 


Corrected Euler 


400 


14322 



Table 1: Computation time in seconds to generate 10^ paths up to T = 1 with = 10 time-steps of the 
fohowing MRC process: k — 1.25/d, c = Id, a = Id, and Xij = 0.7 for i ^ j. 



defined for any parameters satisfying ([3]), the second order scheme given by (H^ and is rather faster 
than the "direct" one presented in Appendix |DJ However, it has a larger bias on our example in Figure ([TJ, 
and their overall efficiency is similar. Nonetheless, both are as expected overtaken by the fast second order 
scheme Let us recall that it is only defined under Assumption which is satisfied by our set of 

parameters. Also, the fast first order scheme given by (j48p requires roughly the same computation time. 

Let us switch now to Figure [1] that illustrates the weak convergence of the different schemes. To be more 
precise, we have plotted the following combinations the moments of order 3 and 1 (i.e. respectively 



E 



kl 



l<i7^i<3 



{X- 



2,3(A^T^)l,3 



(50) 



and E J2i<i^j<d(-^T)i,j ) function of the time-step T/N. These expectations can be calculated exactly 
for the MRC process thanks to Proposition^ and the exact value is reported in both graphics. As expected, 
we observe a quadratic convergence for the second order schemes, and a linear convergence for the first 
order scheme. In particular, this demonstrates numerically the gain that we get by considering the simple 
change of time tp between the schemes (|48|) and (|49|) . Last, the modified Euler scheme shows a roughly linear 
convergence. It has however a much larger bias and is clearly not competitive. 



4 Financial application of correlation processes 

In this section, we focus on the modeling of the dependence between d risky assets. We will denote by 
Si,. . . , Sf their value at time t, and we set log(S't) = (log(S'j ), . . . , log(S'j''))-^. Basically if we exclude jumps, 
we can assume that the assets follow under a risk neutral probability space (fi, T, P) the following dynamics: 

d\og{Sl) = (^r - i(St).,.^ dt + i^tdBt),, 1 < ^ < d (51) 

where r is the interest rate, {Bt,t > 0) is a d-dimensional standard Brownian motion and (St,t > 0) is an 
adapted (E)-valued process. This process describes the instantaneous covariance of the stocks 

{d\og{Sl),d\og{Sl)) ^ {llt)^^Jdt, (52) 

and {^t)i.i is usually called the volatility of the asset SI . 

Modelling directly the whole covariance process (St,t > 0) is not an easy task. This path has recently 
been explored by Gourieroux and Sufana [TH]. Their model has been enhanced by Da Fonseca et al. [S]. 
They assume that (St,t > 0) is a Wishart process and consider a dynamics that is a natural extension of 
the famous Heston model [TB] to d stocks. Since Wishart processes are closely connected to MRC processes, 
we will discuss in detail their model in Section l4.ll In particular, we explain why, in our opinion, it can be 
hardly used in practice when the number of assets is large (say c? > 5). 
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2nd "direct" 




E5:act value 



Figure 1: d = 3, same parameters as for Table [TJ In the left (resp. right) side is plotted (|50p (resp. 



E 



J2i<i^j<d(^' 



in function of the time step The width of each point represents the 95% 



confidence interval (10^ scenarios for the modified Euler scheme and 10* for the others). 



Instead, the common practice of the market is to start with modeling the volatility ctj of each stock SI- 
Then, the dependence is modeled by a correlation process {Ct,t > 0), so that the covariance process is 
defined by 

Sf = diag{al, at)Ctdiag{a] , ...,af). 
Thus, if we consider the following dynamics for the assets 

d\og{Sl) = (^r - i(aj)2^ dt + ali^tdBt),, 1 < ^ < d, (53) 

we get back the same instantaneous covariance given as ()52|) . This bottom- up approach has many advantages. 
Indeed, the modelling of individual stocks is well documented and handled every day by the financial desks. 
Besides, the choice of the volatility model is free and can be chosen at one's convenience. For example, we can 
take a local volatility model (crj = a{t, SD) or a stochastic volatility model such as Heston model, or a local 
stochastic volatility model |3] . The calibration of the single stocks can be thus performed separately, before 
the calibration of the correlation process. There is still few literature that brings on fitting the correlation 
to market data. Today, the only liquid and quoted derivatives that bring on the dependence between assets 
are index options. Recently, Langnau [22] and Reghai [25^ focused on the calibration of some particular local 
correlation models (i.e. Ct = C{t, , . . . , Sf)) to index option prices. Under a slightly different setting, 
Jourdain and Sbai |18_ have also considered this issue. However, there is still up to our knowledge very few 
studies on stochastic correlation modelling. We will discuss in Section 14.21 how MRC processes and possible 
extensions could be used for that purpose. 

4.1 The Wishart stochastic covariance model 

Da Fonseca et al. [9] assume that (St,i > 0) follows a Wishart process 

dEt = {aaa^ + hT.t + T.th^) dt + ^tdWta + dW^ ^/Y^t (54) 
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with a > d — 1, a, 6 G A^c;(R). The vector Brownian motion B and the matrix Brownian motion W may be 
correlated as follows: 

dBt^dWtp+^l-\\p\\ldBt,peR'^ s.t. \\p\\2 < 1, 

where {Bt,t > 0) is a d-dimensional Brownian motion independent from W. With this choice, the couple 
(log(5't), Et) is an affine process, and its characteristic function can be obtained by solving Riccati equations 
(see Da Fonseca et al. [9]). 

From (|54p . we get that there are Brownian motions (31, i ~ 1, . . . ,d such that the volatility square of SI 
follows: 

d 

d{T,t)i,i = {a{aa^)i^i + 2y^^bij{T,t)i,jdt + 2^/ {aa'^)i,i^/ {Y.t)i,id(3l. 
i=i 

The non diagonal elements (St)i.j for i ^ j can be seen as factors that drive the volatility of each stock. In 
its full form, the dynamics of one stock and its volatility is not autonomous. This means in practice that 
it is not possible to calibrate this model separately to single stock market (mainly, European options on 
single stocks). This calibration has to be done at the same time for all the stocks, which is a priori a very 
challenging task: many parameters and a lot of data are involved when the number of assets d gets large. 

Then, one may want to recover autonomous dynamics for each stock in order to calibrate them separately. 
Within this model, the only possible choice is to assume that & is a diagonal matrix. In this case, {^t)i,i 
follows a CIR diffusion and the couple {{l^t)i.i, SI) follows an Heston model [16]: there are independent 
Brownian motions /3* , 7* such that 

dlog(^,') = (r-i^)di+y(^(pV/3j + yr~Wrf7t) 
d{T,t)i^i = (a(a^a)i,i + 2bi^i{'St)i^i) dt + 2yJ {a^a)i^i\J {T,t)i.idPl, 

^^^^ VeLi 4, ' 

In practice, each individual stock could be then calibrated like in the Heston model. Unfortunately, there are 
further restrictions implied by this model and especially that a> d — 1, which is the condition that ensures 
the existence of the Wishart process. This is unlikely because when calibrating Heston to market data, it is 
typical to get values of a around or below 1. When modelling many stocks together (say d > 5), it is then 
not possible to fit conveniently market data because of this restriction on a. 

Last, one of the main feature of the model (|54l) is the Affine property. It allows to obtain the characteristic 
function of the stocks by solving Riccati differential equations. Then, the pricing of European style options 
can be made by using Fourier inversion. This approach is known to be very efficient in a one-dimensional 
framework (see Carr and Madan |6]) and has been used successfully by Da Fonseca et al.[9] to price Best-of 
options with d = 2 assets. However, when the number of assets d is much larger, the Fourier inversion 
requires an integration in dimension d and can no longer be computed quickly. Unless the payoff has a very 
particular structure to reduce the dimension, the pricing by Fourier inversion is no longer competitive with 
respect to Monte-Carlo methods and the Afhne property does not really give a crucial advantage in terms 
of computational methods. 

For all these reasons, we believe that this model can be used successfully in practice for a small number 
of stocks {d = 2, 3) but is instead inadequate to model a large basket. 

4.2 Towards a stochastic correlation model 

Now, we would like to discuss the application of the MRC processes under the framework ([55)1 . The first 
natural idea would be simply to take {Ct)t>n ^ MRCd{Co, k, c, a). With this choice, we would get analytical 
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formulas for correlation swaps. Indeed, a correlation swap between assets SI and {i ^ j) with maturity T 
has, up to a discount factor, the following price: 



dt 



1 



(55) 



To be precise, the payoff of a correlation swap is defined as the average over the period [0, T] of the daily 

and is here approximated by y {Ct)ijdt. Unfortu- 
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nately, up to now, correlation swaps are not quoted on the markets and are only dealt over the counter. It 
is then not possible to get data on their prices in order to calibrate k, c and Co, which would have been very 
tractable thanks to formula ([55]) . 

The only quoted options that bring on the dependence between assets are index options. We have 
been kindly given by Julien Guyon at Societe Generale market data on the DAX index at the 4th October 
2010. He provides us with data on European index option prices as well as parametrized local volatility 
functions a^{t,x) that are already calibrated to options price on each asset. We thus assume in the sequel 
the dynamics ([55)1 for the stocks with ctj — a'^{t, SI). We will also denote 

d 



the index value and will assume constant weights ai such that X)i=i ~ 1- These weights are given in 
Table [2] for the DAX. 



SIEMENS 


BASF 


BAYER 


EON 


DAIMLER 


ALLIANZ 


SAP 


D. TELEKOM 


D. BANK 


RWE 


9.91 


8.03 


7.97 


7.67 


7.35 


6.97 


6.01 


5.6 


5.06 


3.86 


M-RUECK 


LINDE 


BA-IW 


VW 


D. POST 


ADIDAS 


D. BOERSE 


FRESBNIUS-MC 


THYSS. KRUPP 


MAN 


3.23 


3.06 


2.92 


2.26 


2.05 


1.76 


1.65 


1.63 


1.52 


1.47 


HENKEL 


SDF 


LUFTHANSA 


METRO 


INFINEON 


HEIDEL. 


FRESENIUS 


BEIERSDORF 


COMMERZBANK 


MERCK 


1.29 


1.17 


1.16 


1.12 


1.02 


0.94 


0.9 


0.86 


0.84 


0.74 



Table 2: DAX index composition, in percentage (4 October 2010) 



To calibrate MRC processes to index options, it is desirable to reduce the number of parameters. We will 
assume in the sequel with a slight abuse of notation that k = Kid, a = aid with k, a € M.'!^, = li=j + p^i^j 
for p e [~l/{d — 1), 1] and Co — c. Condition ([3]) is then simply equivalent to 2k(1 — p) > {d — 2)a^, but we 
will assume in addition that 

2k(1 - p)>id- l)a^ 

in order to take advantage of the 0(d"^) discretization scheme given by (P^ . Thus, Ct follows the following 
dynamics 

Ct = c + 2K [ {Cs- c)ds + a V / (^Cs - Cse^CsdWsc'ij + e'^dW^ ^Cs - Cse^Cs) , (56) 

and we assume that the Brownian motion W is independent from the Brownian motion B that drives the 
stocks ([5^ . We would like to calibrate such a process to European index option prices. To do so, we first 
calculate the value p such that the constant correlation model li^j + p'^i^j fits the at the money implied 
volatility. Then, we use this value and look at the impact of k and a. Figure [2] illustrates these results. As 
one could expected, this model gives a smile which is not enough sloping and is unable to fit the volatility 
skew. The parameters k and a have no impact on the slope of the smile. 
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Figure 2: In both graphics, the hne with diamonds indicate the market 1 year imphed volatihty. We plot the imphed 

given by model (|56|l with p = 0.67968 and use the parametrization a = ^rf(l ~ p)) with e £ [0, 1]. Left: k — 100 

and e G {0,0.6, 1}. Right: e — 0.5 and k £ {1,20, 100} . These curves are obtained with 10^ Monte-Carlo samples 
and discretization time step 0.025. 

The reason of this too flat smile is that the correlation is not related to the index price, while the market 
would expect that the correlation is high (resp. low) when the index is low (resp. high). There are at least 
two ways to correct this. First, one may assume some dependence between the Brownian motions W and 
_B, which would be very analogous to what is made for stochastic volatility models. This requires to flnd 
an adequate way to correlate these Brownian motion from a financial point of view that is also tractable for 
simulation purposes. We have left this for further research. 

The second way to make depend the correlation in function of the index value is simply to assume a local 
correlation model C{t, It) with C : M+ x — > (^^(R) and keep W and B independent. This approach has 
been considered by Reghai 25 . In Ahdida [I] , it is shown that the following parametric form 

C(t,It)i,3 = li=j +p(i,/t)lj#j, with p{t,It) = max ( — } .7 ,Pmin ) , (57) 

with 7/, 7, Pmin > is vcry tractable to calibrate the index smile at a given maturity. Indeed, it is shown 
that f] tunes the implied volatility at the money, 7 tunes the skew (i.e. the slope at the money of the implied 
volatility), and Pmin tunes the right tail of the smile. The left hand side of Figure [3] shows the index smile 
data and the implied volatility given by this model. Also, an extension of this model with time-dependent 
parameters can be used to fit the index smile for different maturities. 

We want to illustrate now how dynamics such as MRC processes could be used to extend local correlation 
models and add a new source of randomness. Namely, we consider the following dynamics 

Ct = Co + 2k [ {C,-C{s,Q)ds+ [ a{s,Q V (VCs - Cse^CsdWse'^ + e^dW^^C, - C.e^Cs) , (58) 
Jo Jo ^ ^ 

where k > 0, a : x and C{t,I) is defined by ([57]). The limit case a{t,x) = 0, k = +00 

corresponds to the local correlation model. Following the same lines as in the proof of Theorem [T^l we could 
show under some rather mild assumptions on p{t,x) and a^{t,x) that the whole SDE on {Ct, St , ■ ■ ■ , Sf) 
has a weak solution if 2k(1 — p{t,x)) > {d — 2)a^(t, x) for all i > 0,a; > 0, and a unique strong solution if 
2k(1 — p{t,x)) > da^{t,x). To simulate such an SDE, we will simply use the Euler-Maruyama scheme for 
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Figure 3: Left: diamonds indicate market 1 year implied volatility for the DAX. The curve fitting this data is the 
implied volatility given by the local correlation model (|57p with 7 — 8.672568, ri — 0.5007f4 and pmin ~ 0.1. The 
other curve is the implied volatility given by a constant correlation model {Ct)i,j = li=j + plijtj that fits the at the 
money implied volatility: p = 0.67968. Right: With the same values for 7, and pmin, we have plotted for k = fOO 
and £ £ {0,0.3,0.9} the implied volatility given by model (|58p . These curves are obtained with 10^ Monte-Carlo 
samples and discretization time step 0.025. 



the stocks and use our scheme for the MRC process with coefficients k, C(if and a(if on the 
time-step [tf ,tf_^i]. More precisely, we will assume moreover that 2k(1 — p{t,x)) > {d— l)a'^{t,x) and even 
set 

a{t, x) = <^-^^{l- p{t,x j), 

where e e [0, 1] is a free parameter. This choice allows to use the 0{d^) discretization scheme p9)) for the 
MRC process. Starting from the calibrated local correlation model, we have plotted in the right-hand side 
of Figure |3] the effect of the volatility on the index smile. We have chosen a large value for k so that the 
model (f58|) fits the data for e = 0. We see that the volatility of the correlation tends to reduce prices of 
call option on the index. The same monotonicity was already observed in Figure (3] This indicates some 
concavity of the index option prices with respect to the correlation. 

Last, a natural question is to wonder if this is really necessary to sample a whole correlation process. For 
example, we could consider the following onc-dimensional model 



Ct = h=j + Ptli^j, with pq + K {ps- p{s, Is))ds + a - Ps)dWs, (59) 

Jo Jo 

with G [0, 1], K > 0, a > 0. This dynamics would have rather close qualitative features to ([55)) and is 
much less demanding in terms of computational effort. To be fair, as far as index modeling is concerned it 
may be sufhcient to parametrize the correlation matrix by a single parameter pt. The heuristic reason is 
that index options do not really depend on the individual pairwise correlations {Ct)ij but rather depend 
on an average correlation in the basket. Instead, if the aim is then to price and hedge exotic products on 
the dependence, it may be relevant to model all the pairwise correlations. To give a caricatural example, an 
option on the difference of two correlation swaps that pays Jo {Ct)i,jdt— ^ Jq {Ct)k,idt)^ is almost surely 
equal to zero in model (|59p or (|57p . which may basically give an arbitrage. It has instead a non trivial price 
if {Ct)t>o ~ MRCd{Co, K, c, a) or in model 
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Wc' have tested MRC type dynamics on index options because they are the only hquid quoted options 
that bring on dependence. The tractabihty offered by these processes is not really exploited for such options. 
Unfortunately today, there is no quoted options that could give the market view on pairwise correlations. 
However, if an investor has some personal views on correlations between some companies or some industry 
sectors, processes such as MRC can be a relevant tool to take into account these views and price exotic 
products. Generally speaking, modelling precisely the dependence between the stocks in order to get a 
model that prices consistently single-name and basket products is an important challenge in finance, and we 
hope that processes such as MRC may be tool to achieve it. 
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A Some results on correlation matrices 

A.l Linear ODEs on correlation matrices 

Let b e 5£i(M) and k G A^(j(M). In this section, we consider the foUowing linear ODE 

x'{t) = b- {Kx{t) + x{t)K^), x(0) = xe €d{M.), (60) 
and we are interested in necessary and sufficient conditions on k and b such tliat 

Va; e £d(K),V< > 0,x(i) e (61) 
Let us first look at necessary conditions. We have for I < i,j < d: 



c'ijit) = bij - ^ Ki^kXk,j(t) + Xi.k{t)K 



k=l 



In particular, we necessarily have a; - ^{t) — 0. This gives for t ^ 0, I ^ i and a::(0) = Id + p(e^' + ef) that 
bi^i — 2Ki^i — 2pKi^i = for any p e [—1, 1]. It conies out that; 

Kj,/ ^ if I ^ i, bi^i = 2Kj,i. 

Thus, the matrix k is diagonal and we denote ki ~ Kij. We get ^{t) — bij — {ni + Kj)xij(t) for i ^ j. If 
Ki + Kj = 0, we have Xij(t) — Xij + bijt, which implies that bij = 0. Otherwise, Ki + kj ^ and we get: 

= exp {—{ni + Kj)t) H ^-^ — [1 — exp {—{Ki + Kj)t)] . 

Ki-i ~\~ f^j 

Once again, this implies that Ki + Kj > since the initial value x G £d(K) is arbitrary. We set for 1 < i, j < d, 



^3 

" ^ "-J 



if Ki + Kj > 



c.,, = 1, and for i ^ j, c.,j = <^ '''+'^^ ' (62) 

[OifK, +Kj=0. 

We have & = kc + ck and for x = Id, c ^ limt_>._|_oo x{t) G £d(K), and deduce the following result. 

Proposition 19 — Let b e 5(i(M) and k e Ald(M). // the linear ODE ((60)) satisfies (|61l) . then we have 
necessarily: 

3c e £(i(K), 3/ti, . . . , e M, V« 7^ J, Ki + Kj > 0, /t = diag{Ki, . . . , k^) and b — kc + ck. (63) 

Conversely, let us assume that ([M)) holds and b e 5j'(M). We get that Ki — 6^.^/2 > and for t > 0, 
exp{Kt)x{t) exp{Kt) = x + exp{Ks)b exp{Ks)ds is clearly positive semidefinite. Therefore, (j6ip holds. We 
get the following result. 

Proposition 20 — Let ki, . . . , > 0, k = diag{Ki, . . . , Kd) and c € £rf(R). // kc + ck <E '5j'(M) or d = 2, 
the ODE 

x'{t) = k(c - x) + (c - x)k, a;(0) = x G ^^(M) (64) 

satisfies (pT|) . 

Let us note here that the parametrization of the ODE (|64p is redundant when d = 2, and we can assume 
without loss of generality that ki — K2 for which kc + ck £ 5j"(M) is clearly satisfied. 
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Remark 21 — The condition given by Proposition is necessary hut not sufficient, and the condition 
given by Proposition \20[ is sufficient but not necessary. Let c? = 3 and c = I3. We can check that for 
K = (1, |, —\), holds but (pT|) is not true. Also, we can check that for k — (1, 1, —\), (pT|) holds. 



Lemma 22 — LetK^,K^ be diagonal matrices and c^ , c^ G such that k^c^+c^k^+k^c^+c^k'^ ^ 

Then, the ODE 

x' = n^{c^ - x) + (c^ - x)k^ + K^{c^ - x) + (c^ - x)k^ 
satisfies (j6ip . Besides, x' — k{c ~ x) + {c — x)k with k = + k'^ £ S'^(R) and c G defined by: 



1, and for i ^ j, ct. 



if Ki + Kj — 0. 



Proof : Since b = k^c^ + c^k^ + + (?k} E iSj"(M), (pT|) holds for x' — b — kx + xk. Then, we know 
by that c is a correlation matrix. □ 



A. 2 Some algebraic results on correlation matrices 

Lemma 23 — Let c £ and 1 < i < d. Then we have: c — ce^^c £ S^{M.), (c — ce\c)i^j ~ for 

l<3 <d, (c - cejjc)''' = cW - c\dY and: 



Besides, if c £ - d{d Y £ 5+4 (R). 



= 0. 

3 



Proof : Up to a permutation, it is sufficient to prove the result for 1 = 1. We have 

c-ce\c = ^[1] ) = aca^, with a = "Jt! ) ^ ^^^^'^^ 

Besides, we have Rk(aca^) = Y^]^{a^/c) = d—l when c G which gives — c*(c')^ G S'^^^{ 

Lemma 24 — Let c £ (^^(R) a'^'^ I <n < d. Then — ,/ce^^^/c £ 5j"(K) and is such that 



Id - \/ce^'Vc ^ Id - \fce^'sfc. 



Proof : The matrix (%/ce2Vc)i,j = (Vc)i,n(\/c)j\„ is of rank 1 and I]^=i(Vced Vc)ij (a/c)^/™ = (Vc)i,ri since 

n = = 1. Therefore ((■\/c)i,n)i<i<rf is an eigenvector, and the eigenvalues of 1^ — ^/ce^^/c are 
and 1 (with multiplicity d—l). □ 



Lemma 25 — Let q £ iSj"(M) be a matrix with rank r. Then there is a permutation matrix p, an invertible 
lower triangular matrix nir £ QriJ^) and k^ £ AAd-r^riJ^) such that: 



T T 
pqp — mm , m 



mr 
kr 



The triplet {mr,kr,p) is called an extended Cholesky decomposition of q. 
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The proof of this result and a numerical procedure to get such a decomposition can be found in Golub 
and Van Loan {[13], Algorithm 4.2.4). When r ^ d, we can take p — Id, and is the usual Cholesky 
decomposition. 

Lemma 26 — Let c € <Ld{^), r — Kk{{cij)2<ij<d) o,nd {mr,kr,p) an extended Cholesky decomposition 
of {ci.j)2<i,j<d- Wesetp^(^^ ^j. 

where c\ E W , with {ci)i — ijp" cp)i^i+i for 1 < i <r. We have: 

c = pmcm'^p^ and c G €d (K) ■ 

Proof : By straightforward block-matrix calculations, on has to check that the vector c^''' G 
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rrir 







{0 












{m-'cir 





and c — 




Ir 







K 





Id-r-1 



pd-(r+l) 



defined by (cj^' )i = {p cp)i_i for r + 1 < i < d is equal to krm,. c\. To get this, we introduce the matrix 
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^0 




Id-r-1 













and have q ^ = 
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—krm~^ 


Id~r-1 



Since the matrix 
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{m-'c^r 


/ r.d 


— krnij. ^Ci)'^ 




Ir 







'^^d i. —1 r 











-IT / -1\T 

q p cp{q ) 



is positive semidefinite, we have c^' — krm~^c[, -i r ^ i'' 

'Tij. Cy If 



e 5++i(K) and thus c e ^^(IR)- □ 



B Some auxiliary results 

B.l Calculation of quadratic variations 

Lemma 27 — Let {J-t)t>o denote the filtration generated by (Wt,t > 0). We consider a process (Yt)t>Q 
valued in Sd{R) such that 

d 

dYt = Btdt + J2iA7dWte-d + eyw;'{A^f), 

n=l 

where (A")t>o, (i?t)t>o continuous (J-'t) -adapted processes respectively valued in A^rf(M), and iSd(M). 
Then, we have for 1 < i,j, k, I < d: 

d{Y,,„Yk,i)t = [l^=k{Al{Alf),,i + l,=i{Al{Alf),,k + l,=k{Al{AlfU + l,=i{Al{Alf),,k] dt (65) 
Proof : Since {A^dWtc''^),,, = lj=„(A^ dW^t)^J and (e2rfWt^(A?)^),,, = l^=n{A\dWt),,^, we get: 

d 

d{Yt)^,, = {Bt)^^,dt +Y.(^tkn{dWtUj + {Al)J,,,{dWt)r^,^■ 
n=l 

Then, d{Y,,,,Yk,i)t = [lj=i E'=i(^D.,4^D/c.n + 1,=^ E'=i(^D.,«(A^');,„ + 1.=; E'=i(^j).^n(Aj)fe,„ 

+ li=fe E«=i(^t)j,«(^t)i,n di, which precisely gives ([SS]). □ 
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Lemma 28 — Let us consider x G (^^(R), and {Xt)t>a a solution of the SDE Let t denote the stopping 
time defined as t — {t>Q^Xt ^ £^(R)}. Then, there exists a real Brownian motion {Pt)t>o such that for 

0<t<T, 

d(det(Xi)) 



- Tt[X^Hkc + ck- id-2)a'^)]dt~Tri2K + a^)dt + 2jTr[a^{Xf-^ - Id)]dl3t, (66) 
d\og{det{Xt)) = TT[Xf\KC + ck- da^)]dt ~ Tr(2K - a^)dt + 2 Jxr {a^{X^^ - Ld)] df3f (67) 



Proof : First, let us recall that yi,j,k,l G {l,...,d}, Wx £ 5j''*(M) dij det{x) = (adj(a;))ij- = det{x)x^j , 
dk,idi J {det{x)) = det{x){xY^x~j—x'[jX~l). Since x is symmetric, we have in particular that 3fc.;9ij (det(x)) = 
if i = / or j = fc. Ito's Formula gives for t < t: 



= ixr\AXtkj + l J2 {{X^\AXr\4-iXt-\k{X^')u){diXtk,,diXt)k, 



d(det(Xt)) 

l<k,l<d 

On the one hand we have 

d 



{X^%,jd{Xtkj = Tl[X^^{KC + CK)]dt-T^:{2K)dt + 2Ya^^^ 

l<i,j<d i=l 



X^'el^dWjJXt - Xtel^Xt 



On the other hand we get by ©: 

J2 {{Xi\,{X^\,i - {Xr\k{Xi\i) {d{Xtk,,d{Xt)kd) 



l<i,j<d 
l<k,l<d 



J2 {iXr\AXr')k,i - iX-\kiX^'hi) X {a'^l,=kiXt - Xte^XtU 



l<i,j<d 
l<k,l<d 

+apj=i{Xt - Xte\Xt\k + a?l.=i(^t - Xte\Xt)^^u +a^l,=fe(Xt - Xte\Xt)^^i) 

= E f E - Xte'dXtkk {{xi%^,{xi^k, - (xr^).,fc(xr^),j) 

j=l \\<:i,k<d 

+E f E «?(^* - ^te^^tk; ((xr^).,,(xr^)M - (xr^)M(^r'k/) 

d 

= 2 ^ [{X, - X,e:,Xt)Xi^e:,X^'] - {X^%,^Ti [{X^ ~ X,e\Xt)Xl^\ ) . 

Since e ^(M), we obtain that Tr {(Xt - Xtej;Xt)Xt-iej;Xt-i] = {X^%^^~\ and Tr [X^-^Xt - Xte\Xt 
d — (^t)i,i = d — 1. We finally get: 



d{d&X{Xt)) 
det(Xt) 



{d-2)a^)]dt~Tr{2K + a^)dt + 2Y 
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Now, we compute the quadratic variation of det{Xt) by using (jS]): 
d{detiX))t 



det(Xt)2 



l<i,j<d 
l<k,l<a 



Wili=i{Xt - Xte'^Xt)j.k + a\\,^k{Xt - Xte\Xt),,i\ dt 

d 

^Y^afTr [X^^e^X^HXt - Xt^^X^)] dt 



i=l 
d 



= 4^a,2((X("^)j,, - l)dt = 4[Tr(a^Xt"^) - TT{a^)]dt. 

i=l 

It is indeed nonnegative: we can show by diagonahzing and using the convexity of a; i--> 1/a; that x^l > 
^/xi^i = 1. Then, there is a Brownian motion > 0) such that (|66|) holds (see Theorem 3.4.2 in [H]). □ 



Proposition 29 — Let k,6,ri > 0. For a given x G [—1.1], let us consider a process {X^)t>o, starting from 
X, and defined as the solution of the following SDE 

dX^ ^k{9~ )dt + r] v/1 - {XfYdBt , (69) 

where {Bt)t>o is a real Brownian motion. Then there exists a positive constant if > 0, such that 

yt > o,yx e [-1, 1], E [{x^ ~ x)^] < Kt^ 

Proof : For a given x £ [—1, 1], we set /^(y) = {y ~ x)'^. If we denote L the infinitesimal operator of the 
process Xf, then we notice that f^{x) ~ Lf^{x) = 0. Besides, {x,y) G [—1,1]^ i— ^ i^/^(j/) is continuous 
and therefore bounded: 3^ > g, Va;, y G [-1, 1], |LV(y)l < 2i^. (70) 

Since the process {Xf)t>Q is defined on [—1, 1], we get by applying twice Ito's formula: 



nrixt)] 



■t rs 



E [L^f^X^)] duds. 



From ([70)) . one can deduce that E [L^/^(X|^)] duds < Kt^, and obtain the final result. 



B.2 Some basic results on squared Bessel processes 

Lemma 30 — Let (3 > 2 and Zt — z + /3t + 2 \fZldBs be a squared Bessel process of dimension j3 starting 
from z > 0. Then we have 

/■* ds f^°° ds 

P(Vt > 0, / — < cx)) = 1 and / — = +00 a.s. 
Jo Zs Jo Zs 



Proof : The first claim is obvious, since the square Bessel process does never touch zero under the con- 
dition of /? > 2. (see for instance [17], part 6.1.3). By using a comparison theorem {\/t > 0,Zt < Z[ 
a.s. if /3 < /3'), it is sufficient to prove the second claim for /3 G N. In this case, it is well known that 
{W} + \fz)^ + 'Y%,=2^t)'^ follows a square Bessel process of dimension n, where (W^/'',i > 0) are indepen- 
dent Brownian motion. By the law of the iterated logarithm, lim sup(_^_|_g^ 2tiog(iog(t)) ~ which gives the 
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desired result since 



dt 



tlog(log(t)) 



Lemma 31 — Let/3>6. Let Zt ~ 1 + l3t + 2 j^^/ZldBs he a squared Bessel process of dimension 13 starting 
from 1 and (j){t) = -^ds. Then we have 

nnt)] = t + ^^t" + o(t'), E[0(<)2] = t' + oit^), E[cj)itf] = o(t3). 



Proof : For a fixed time t > 0, the density of Zt is given by: 

z > 0,p{t,z) = 

k=0 ■"■ ^"V'>'-r2 



1 



fc! 2tT{k+^)2t 
Let us consider tliat 7 G {1,2,3}, then all negative moments can be written as 

+00 



E 



1 



1 



We have 



fc=0 

1 



{k+§-l) k+l 2(fc+2)(fc+l) 

+ 00 



E 



1 



E 

fc=0 



r(A^+§) ^' (fc + f-l)x...x(A: + f-7) 
L O(-p-), which yields to the following expansion: 

2 o-TT^ 1 \k+3\ 



+00 



(fc+1) 



fc=0 



(fc + 2) 



= 1 - iP - A)t + Oit^) 



(71) 



The first equality is thus obtained. We use the same argument to get: 

+00 „-A 1 Nfe+2 / +00 



E 



E 



1 



Zf 



E 



(fc + 2) 



fe=0 



(fc + 3) 



= 1 + 0(t) 



\fc=0 



(fc + 3)! 



0(1). 



(72) 



By Jensen's inequality, one can deduce that E 



(/o 



* ds_ 



< t'^E 



Jo 



PIP 



. Thanks to the moment expansion 



in (|72p . we find the third equality. Finally, by Jensen's equality, we obtain that 



ds 



< tE 



1 ds 









■ /■* ds ' 


= tE 


Uo {ZsY\ 


- 2tE 


[Jo iZs)\ 


= 


- 2t^ + ^ 


-o{t')-- 





+ t' 



It yields that 
E 





■ 1 ■ 








■ 1 






■ 1 " 


[a* 




ds^ 


= E 


[(/ 


1 


ds^ 


-t^ + 2t [ E 














Zs 




Jo 


.Z's. 



ds^t^ + 0{t^). 
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C A direct proof of Theorem [6] 

Proof : From ^ we have 2Li = -aLf + Lf, with: 

l<j<d l<j,k<d 

We want to show that LiLj — LjLi for i ^ j. Up to a permutation of the coordinates, Li and Lj are 
the same operators as Li and L2. It is therefore sufBcient to check that L1L2 = L2L1. Since L1L2 = 
Lf^ — a(Lf + Lf^ L2) + a^Li L2 ^ it is sufficient to check that the three terms remain unchanged 
when we exchange indices 1 and 2. To do so we write: 

Ll = Y ('^{'J} ^ + 2 XI (^{2,0 - X{l^2}X{l,^})^{la}^{l^^} + (1 - x\^ 2})d{l,2} 

3<i,j<d 3<i<d 

^2'= Yj (^^{fe,/} - a;{2,fc}a;{2,;})9{2,fc}9{2,/} + 2 ^ {x{i,i} - X{i^2}X{2d})d{i^2}d{2j} + {I - xl^ 2})d{i.2} 

3<kd<d 3<l<d 

if = ^{1,2}'9{1,2} + Y if = 2:{1,2}C'{1,2} + X! 2;{2,;}l9{2,Z}- 

3<j<d 3<;<d 

By a straightforward but tedious calculation, we get : 

Lf^Lj^ = X (^{ij} - 2:^{i,i}2^{ij})(^{/c,;} - 2;{2,fe}a;{2,i})9{i_i}9{i,j}9{2,fe}i9{2,/} 

3^i,j,k,l<d 

^ ^ ^ 

T 

+ - 2;{i,,}a:{i,j}) (2a{i^2}9{2,.}a{ij} + 2a{i,2}a{2,j}5{i,,}) 

2 

+2 Y - xii^qxiijy){xiijy - X[i^2}X{ij})dii^2}d{2d}d{i,i}9{i,j} 

3<i,jd<d 

" V ' 

3 

+ Y (^{ij} "^{l,i}^{l.j})(l~^{l,2})^{l,i}^{lj}^{l,2} 

3<i j<d 

4 

+2 ^ (2;{2,i} - X{ia}X{i.i}){xik^q - X{2,k}X{2d})d{2,k}d{2,l}d{i^2}d{l,t} 
3<i,k,l<d 



( 



+4 Y (2;{2,i}-a:^{i,2}a;{i,i}) 



3<i<(i 



^{l,2}C'{2,i} - XI 2;{2,/}9{l,2}9{2,;}a{i,,} + X -a;{i,2}2;{2,/})9{i^2}5{2,i}^{l 

/ 



3</<c; z<i<d 



"2 X! (^{2.*} ^ 2;{i,2}a;{i,,}) 



3<i<d 



-2a;{i,2}5{i 2}5{i.«} + (1 - a^{i,2})5{i,2}^{i,i} 

V ' " V 

8 9 



- Y (^{fc,/} - 2;{2,fc}a;{2.m})(l - a;{i_2})5{2,fe}9{2,m}a{i^2} 
3<k.l<d 
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+(l-a^{i.2}) 



{2,1} 



3<l<d 



V 



10 



+ (1 - a;{l,2})^{l,2}((l - X{i^2})d{l,2}) 



In this formula, the terms n are aheady symmetric by exchanging 1 and 2. The terms n are paired with 
the corresponding symmetric term. To analyse the terms h, we have to do further calculations. On the one 

hand, 



10 + 5 = ^ 4x{i,2}(a;{i,2}a;{2,o - 2;{i,i})^{i,2}^{2,o 
3<;<<i 



3<l<d 



are symmetric together. On the other hand we have 

2 + 6= Y {"^Hhj} - "^^{hiy^hj} - 'ix{2,i}X{2,j} + 4x{i,i}X{2,j}a;{i,2}} ^^l^2}^{l,i}^{2J}, 

l<i,j<d 
i7^1,2,j>^l,2 

which is symmetric. 

Now we focus on L^L^ + Lf^Lf . We number the terms with the same rule as above, and get: 

L^L^+L^L2= Y/ ~ X{2,l}X{2,k})X{l,2}d{2,l}d{2,k}d{l,2} 

3<k,l<d 

V ' 

1 

+ 2 ^ X[i^2}iX{l,l} - X[i^2}X{2,l})dl^,2}d{2,l} - 2 ^ X{i^2}X{2,l}d[2,l}d{i^2} + X{i,2}d{i,2} [{1 - xli^2})9ll,2}} 



3<l<d 



3<l<d 



+ Y -^{l.^lG^f;.*:} - ■'^{2.fr}-'^{2.;})9{2.fr}9|2./}9{l,,} + 2 Y •'J^{1.»}C'^{1-;} - ■^{1.2}-'^{2./})9{1.2}9|2./}9{i.,} 



3<t./,-./j^(/ 



2 5Z a;{i,i}5{i,2}5{2,i} + Y " ^{l,2})^{l,2}5{l,i} + Y (^{ij} " ^{l,i}^{l,j})^{l-2}^{l,2}5{l,i}5{l,j} 



3<i<d 



3<i<c! 
' ^ 



3<i,i<d 



7 8 1 

Y ^{2,«}(^{i,j} - a;{i,i}a;{ij})a{i,j}a{ij}a{2,i} + 2 ^ a;{i,2}(a;{2,i} - X{x,2-)X{x,i-))d\x,2}d{\,i) 

3<i,j,l<d 



3<i<d 
' ^ 



+ 2 X] '^^{2,i} -a;{l,i}X{i,2})5{l,i}5{l,2} +2 Y X{2,l}{X{2,i} - X{i^2}X{l,i})d[i^iyd[i^2}d{2,l} 



3<i<d 



3<i,l<d 



(1 - •'■{l.-'})'^''U-'}i ' {l.-"/''{.L2}} + XI (1 ~ 4l>2})^{2,!}^{l,2}%,;}- 



3<l<d 



Therefore, Lf^Lf + L2 is symmetric when we exchange 1 and 2. Last, it is easy to check that 
LfLf* = L2L1 , which concludes the proof. □ 
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D A direct construction of a second order scheme for MRC pro- 
cesses 

In Section |31 we have presented a second order scheme for Mean- Reverting Correlation processes that 
is obtained from a second order scheme for Wishart processes. In this section, we propose a second order 
scheme that is constructed directly by a splitting of the generator of Mean-Reverting Correlation processes. 
As pointed in (|^ . it is sufficient to construct a potential second order scheme for MRCd{x, ^^e.\, Id, e^; t). 
Thanks to the transformation given by Proposition [9l it is even sufficient to construct such a scheme when 

{x)2<i,j<d = Id-l- 

Consequently, in the rest of this section, we focus on getting a potential second order scheme for 
MRCd{x,^^e\,IdTe\]t), where (a;)2<ij<rf — Id-i- By (1^ . the matrix a; is a correlation matrix if 
Sf=2 ^1 i — 1- Besides, the only non constant elements are on the first row (or the first column) and 
the vector {{Xt)i,i)2...d is thus defined on the unit ball B: 

D= |a;e]R''-\ < l|. (73) 

With a slight abuse of notation, the process {{Xt)i)i...d~i will denote the vector {{Xt)i^i+i)i...d-i- Its 
quadratic covariance is given by d{{Xt)i,d{Xt)j) = (li=j — {Xt)i{Xt)j)dt, and the infinitesimal generator 
of MRCd{x, ^^e]^,Id,e\) can be rewritten on D, as 



2 

1=1 l<i,j"<d-l 

One can prove that the following stochastic differential equation 



d~l ^ 

+9 Yl ~ x^Xj)d,dj. (74) 



VI < j < d - 1, dMl = — H- Mi 



\ 



1 - YiMlYdBl + (1 - {Mlf)dBl+' -Ml Y MldBt' 

3=1 l<3<d-l 



is associated to the martingale problem of L^, where {Bt)t>o denotes a standard Brownian motion in dimen- 
sion d. By Theorem (|12p . there is a unique weak solution {Mt)t>o that is defined on D. 

The scope of this section is to derive a potential second order discretization for the operator L^, by using 
an ad- hoc splitting and the results of Proposition [181 We consider the following splitting 



L^ = C^ + J2 (75) 

m— 1 

where we have, for 1 < m < d — 1: 

\ i=l ) \<l,k<d~\ 

£'"+1 ^ \\ E ^fc5fe + (l-2;m)'C-2a;™(l-x2j xkdud^^ ^ x^xix^d^di 



l<kyim<d-l l<fc^m<d-l l<kjtm<d-l 

l<l^m<d-l 



Thanks to Proposition I18[ it is sufficient to focus on getting potential second-order schemes for the opera- 
tors 
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D.l Potential second order schemes for C? , . . . .C^ 

All the generators I — 1, . . . ,d — 1 have the same solution as up to the permutation of the 

first coordinate and the l-ih. one. It is then sufficient to focus on the first operator . By straightforward 
calculus, we find that the following SDE 

d{Xt)i ^ {I ~ {Xt)l)dBt , ^2<i<d^ld{Xt\^-{Xt)^[^ + {Xt)idB^ , X^,=xe^}, (76) 

is well a solution of the martingale problem for the generator . The SDE that defines {Xt)i is autonomous. 
Since a;i £ [— 1, 1], it has clearly a unique strong valued in [—1, 1]. It yields that the SDE ([75)1 has a unique 
strong solution on . To prove that {Xt)t>a takes values in B we consider Vt — '^i=i{Xt)i ■ By Ito calculus, 
it follows that 

dVt = (1 - Vt){l - iXt)l)dt + 2{Xt)i{l - Vt)dBt. 

Thus, 1 — Vt can be written as a stochastic exponential starting from 1 — Vb > and is therefore nonnegative. 
We now introduce the Ninomiya-Victoir scheme for the SDE (|76p . 

Proposition 32 — Let us consider a; G D. Let Y be sampled according to P{Y = \/3) = P{Y = — VS) = ^, 
so that it fits the first five moments of a standard Gaussian variable. Then Xf = X^{^,X^{^/iY.X^{^Tx))) 
is well defined on D and is a potential second order scheme for the infinitesimal operator , where: 

Vt > 0, e ©, x) = —=^l£==^ V2 < / < d - 1, Xf{t, x) = 



Ve2*x2 + (l-a;2)' " " ' ^ ^ e^^^ xl + {I - xlY 

VyeK,VxeD, X^(y,2;) = ^J|i±^4-7T^' V2<Z<d-l, Xl{y,x)- 



e2y(l + a;i) + (l-a;i)' " " ' ( ^ e'^y{l + xi) + [l - Xi)' 

Proof : The proof is a direct application of the Ninomiya-Victoir 's scheme [241 and we introduce the 
following ODEs: 

dtX^^{t,x) - X^^{t,x){l - {X^{t,x)f), V2 < / < d - 1, dtX^{t,x) - -X°{t,x){X°{t,x)f 
dyXl{y,x) = (1 - (Xi(y,x))2), V2 < / < d - 1, dyXl{y,x) = -Xl{y,x)Xl{y,x). 

These ODEs can be solved explicitly as stated above. We have to check that they are well defined on D. 
This can be checked with the explicit formulas or by observing that 9t(^f~]^(X°(t, x))^) — 2(X°(t, x))^{l — 
EtliXi'it^x))^), dt{T,tl{Xl{t,x))^) = 2Xl{t,x){l-J2f=liXl{t,x)y). Last, Theorem 1.18 in Alfonsi [4] 
ensures that X^ is a potential second order scheme for £'^. □ 



D.2 Potential second order scheme for 

Let {Bt)t>o be a real a Brownian motion. We consider the following SDE; 



VI < i < d - 1, d{Xt)^ = (At), 



1- ^(At)^dBt, Ao = a:eD (77) 



Its infinitesimal generator is and we claim that it has a unique strong solution. To check this, we set 
Zt — yYlt^ii^t)!, - By Ito calculus, we get that the process {Zt)t>o is solution of the following SDE 



dZt = Zt\ll~ZjdBu Zo = 



-1 



\Y.^l (78) 
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Since the SDE ([75)) satisfies tfie Yamada-Watanabe conditions (Proposition 2.13, Cfiapter 5 of [IH]), it has 
a unique strong solution defined on [0, 1]. If Zq = 0, we necessarily have Zt ^ and thus {Xt)i = for any 
t > 0. Otherwise, we have by Ito calculus d\n{{Xt)i) = d\n{Zt), and then 



[^Zt otherwise. 



(79) 



Conversely, we check easily that ([7^ is a strong solution of (|78l) . which proves our claim. The explicit 
solution (l79l) indicates that the SDE ([78]) is one-dimensional up to a basic transformation. Thanks to the 
next proposition, it is sufficient to construct a potential second order scheme for Zt in order to get a potential 
second order scheme for (l78l). 



Proposition 33 — Let us consider x G D, and Z^ denote the second potential order scheme for {Zt)t>Oy 
starting from a given value z G [0, 1]. Then the following scheme Xf 







Vl<i<d-1, (Xf)., 



= otherwise, 



is a second potential order scheme for which is well defined on . 



Proof : For a given a; € E' and / G C°°(D), let (Xf )t>o denote a process defined by ([7^ and starting from 
X G D. It is sufficient to prove that 



E[/(xn]-iE /(In 



The case where cc = is trivial, and we assume thus that ^1 > 0. Let / G C°°(D). Wc define g^ : [0, 1] — >■ 



by Vy G [0,1], g^{y)= f{- 



—====y ... , ^\ »). Since for every 1 < « < d — 1. i , ^ , 

follows we can construct from a good sequence of / a good sequence for g^ that does not depend on x. By 
the defintion of the second potential scheme, there exist positive constants K > Q and ry > 0, depending only 
on a good sequence of / such that G [0, rj] 



< 1, it 



E 



which gives the desired result. 



E 



< Kt^, 



□ 



We now focus on finding a potential second order scheme for {Zt)t>Q- To do so, we try the Ninomiya- 
Victoir's scheme [25 and consider the following ODEs for z G [0, 1], 



Vi> 0, atZo(t,z) = Zo(t,z)(Zo(t,z)- -), VxgK, c>,Zi(x,z) =Zi(x,z)Vl-^i(a;,2)2. 

These ODEs can be solved explicitly. On the one hand, it follows that for every i > and z G [0, 1] 

z exp(— 



Zo(t, z) = 



^l-2z2(l~exp(-i))' 



On the other hand, we get by considering the change of variable ^Jl — Z^ that for every a; G M and z G [0, 1], 

( 22cxp( — a:) 



Zi{x,z) 



, , if a; < 1 lnCi+:^^S^ 

l-y]^l2+oxp(-22:)(l + /T^) - 2 ^'-^\i^^TZ^ h 

1 otherwise. 
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Then, the Ninomiya-Victoir scheme is given by Zo(t/2, Zo(t/2, z))), where y is a random variable 

that matches the five first moments of the standard Gaussian variable. Unfortunately, the composition 
Zo[t/2, Zi{\/tY, Zf){t/2, z))) may not be defined if z is close to 1. To correct this, we proceed like Alfonsi [4] 
for the CIR diffusion. First, we consider Y that has a bounded support so that Zo{t/2, Zi{^/iY, Zo{t/2, z))) 
is well defined when z is far enough from 1 (namely when < z < K{t) < 1 with K{t) = 1 + 0{t)). When 
the initial value z is close to 1, we instead use a moment- matching scheme, and then we prove that the whole 
scheme is potentially of order 2 (Propositions |34] and l35t . 



D.2.1 Ninomiya-Victoir's scheme for {Zt)t>Q away from 1 

Proposition 34 — Let us consider a discrete random variable Y that follows ¥(Y = VS) = F(y = — V3) = 
and P{Y = 0) = |, so that it matches the five first moments of a standard Gaussian. 

• For a given z € [0, 1], the map z ^ Zo{t/2, Zi{VtY, Zo{t/2, z)))) is well defined on [0, 1], if and only if 
z G [0, -fC(t)], where the threshold function K(t) is given in ()81|) . 

• For a given function f e C°°([0,1]), there are constants rj,C > depending only on a good sequence 
of f such that yt e [0,77], Vz e [0,K{t)], 



Zo(t/2, Z^iVtY, Zoit/2, z)))) - ( /(z) - tLzfiz) + yL|/(z) 



< Ct^, 



(80) 



where Lz is the infinitesimal operator associated to the SDE ([78]). 
For every t >Q the function K(t) is valued on [0, 1] such that 



K{t)^ 



2 - e-*/2 



-*/2+2(l-i^(t,V3)2)(l-e-*/2) ^ ^ 



(81) 



with \fy e M+ D{t, y) 



-tji 



(l+e 



e-2VTB + l+ 



l-<!-*/2 



t/2 y- I 



Proof : The main technical thing here is to check the first point. Then, (1501) is a direct consequence of The- 
orem 1.18 in Alfonsi W. By construction, we have 2'o(t/2, z) e [0, 1] <^ z < . We conclude that 

Y^2-oxp(t/2) 

the whole scheme Zo(V2, Zl{^ftY, Zo[t/2, z)))) is weh defined on [0, 1], if and only if Zi{^ftY, Zo{t/2, z))) < 
I ^ By slight abuse of notation, we denote in the following Zo(i/2, z) by the shorthand Zq. Let us 

(t/2) 

assume for a while that we have: 

^1 - Zl{\ + e-2VtV) + g-2v^y - 1 > 0, a.s. (82) 
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It yields then to 

Z^{t/2,Z^{^tY,Zo{tl2,z)))) e [0,1] 



1- Zi{VtY,Zo{tl2,z))) 



> 



'1 -e-*/2 
2-e-*/2 



\ 



> 



e-2v^^(l + y/T^) - (1 - y/T^) 
^e-2^/*v■(l + ^T^) + (1 - yT^)y V 



> 



By in} y/T~^(e-2Vt^ -!) + ! + 6-2^*5- " V 2 - e' 



-t/2 



'l-e-*/2 
2 - e-*/2 

(83) 



1 


— e 


t/2 




— 


t/2 




— 


t/2 


2 


— 


t/2 



■.= D{t,Y). 



2(1 + 



-t/2 ' 



We can check that the mapping D : {t, x) E K+ x R D{t, x) = —1 -\ ,' is non 

e-2vT. + i+^i^£Zi^(l_e-2vT.) 

decreasing on x, and D{t, x) < 1. Since Y G {— \/3, 0, Vs] , it yields thus that the last condition is equivalent 
to: 



Zo(t/2, z)< Jl- D{t, \/3)2 ^ z< 



i-Dit,V3y 



e*/2 + 2(1 - e-*/2)(i _ i:)(i, 73)2) 



(84) 



Conversely, if dHU) is satisfied, we can check that D{t,Y){l + er'^^^) + e-^'^^ - 1 > 0. Therefore (15^ 
holds. To sum up, when z G [0, i^(i)], we both have Zo{t/2, z), Za{t/2, Zi{ViY, Zo{t/2, z)))) e [0,1]. 

Last, it remains to compute the limit of (1 — K{t))/t. First, it is obvious that lim(_j.o K{t) — 1. We can 

check that Jl~D^{t,V3) = ^„ ^ a'"^,"!!,., — = 1 + - ^(1 + VS)) + o{t), and therefore 



1 - 



^/l-D{t,V3'f 



Vl+2(l-r>(t,V3)2)(l-e-'/2) 



\/2-e-'/2+7l-e-'/2(l-e-v^^) 

= t(^(l + V3)) + o(i). It yields that limt_,o ^^^f^ = ^(1 + VS) V i. 



□ 



D.2.2 Potential second order scheme for {Zt)t>o in a neighbourhood of 1 

Let {Zt)t>Q be solution of the SDE ([751) . By ltd calculus, its moments satisfy the following induction: 

Vfc > 2, E [Z^] = (z^ - 1^ M^e-^^E [^,^+2] rf,^ exp(fcllt). 
We obtain first that E [Zf] =7^ + 0{t), then E [Zf] ^ z^ + &zH{l - z'^) + 0{t^) and last 

E [Z2] = z2 + tz\l - z2) + ^lz\l - z2)(l _ 6z2) + 0(^3). (35) 

Moreover, by straightforward calculus, one can check that if t < | and for every z G [0, 1] 

z2+te2(l-z2) + ^z2(l-z2)(l_6^2) < ^ ^ ^2 (1 - z2) + ^ z2 (1 - z2 ) (1 - 6z2) > 0. (86) 

Since E(Zt) = z, the right hand side of (|86p corresponds to the asymptotic variance of Zf . To approximate 
the process {Zt)t>Q near to one, we use a discrete random variable, denoted by Z^, that fits both the 
exact first moment and the asymptotic second given by (|85p . We assume that Z^ takes two possible values 
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< z+ < z~ < 1, with probability p{t,z) and 1 — p{t,z) respectively. We introduce two positive variables 
(rn'^,m^), defined as z+ = z + m"*" and = z — . Since we are looking to match the moment, we get 
the following equations: 



E =z m+p{t, z) = ■m~{l - p{t, z)) 

\ZM = z2 + ^2(1 - z2) + ^^2(1 _ ^2)(i „ g^2) ^ (m+f^ = tz^il - Z^) + t^Z^l - Z^)il - 6z^) 

(87) 

We choose 



E 



z{l — z) and then have p{t, z) = 1 — 



1 



The random variable is well defined on [0, 1] if and only if z+ < 1 and z^ > 0, which is respectively 
equivalent to z(l — z) < (1 — z) and t{l + z)(l + |(1 — Gz^)) < 1. By straightforward calculus, we can check 
that these conditions are satisfied. Since 1 — K{t) = 0{t) by Proposition [Ml we deduce that there is C > 

such that 

e [0, ?l,Vz e \K(t), 1], Vg e N*, E [(l - Z'y] < CH" (88) 
5 L J 



Proposition 35 — Let U ^ U{[0, 1]). The scheme = z+l{(7<p(j 2)} + ^~ l{;7>p(t,2)} 0. potential second 
order scheme on z € [K(t),l]: for any function f G C°°{[0,1]), there are positive constants C and 11 that 
depend on a good sequence of the function /, such that 



Vie [0,77A-], yze[Kit),l] 
5 



E 



-f[z)-tLzf{z)--{Lzff{z) 



(89) 



where Lz is the infinitesimal operator associated to the SDE (j78p . 



Proof : Let us consider a function / e C°°([0, 1]). Since the exact scheme is a potential second order 
scheme (see Alfonsi [4]), there exist then two positive constants r] and C, such that Vi G [0, /i], Vz G 
[0,1], \R[f{Z^)] - f{z) - tLzf{z) - tL{Lzff{z)\ < Ct^. We conclude that it is sufficient to prove that 
Vz G [K{t),l], |E[/(Zf)] -E[/(^f)]| < Ct^, for a constant positive variable C. By a third order Taylor 
expansion of / near to one, we obtain that 



VzG [0,1], /(z)- /(l)-/'(l)(l-z) + 



[l-zf 



/"(I)) <||/^=^^||oo(l-z) 



Thus, there is a constant C > depending on a good sequence of / such that 

|E [/(Zf )] - E[/(Zf )]| < C (e[(1 - Ztf] + E[(l - Zif] + E [(1 - Zff] - 



(1 - Zf)^ 



By IMl), the first term is of order 0{t^). The last term is equal to |E[(Zf )2] - z^ - tz^{l - z^) - (1 - 
z2)(l - 6z2)| and is also of order 0{t^) by (US]). Last, we have by Ito calculus that Vg > 2, E[(l - Zf )«] < 
(1 — zY + (/(q — 1) JJe[(1 — ZlY~^]ds. By induction, we get that there is a constant Rq > 0, such that 
Vz G [Js:(<), 1], E[(l - Z'iY] < Rqti, which finally gives the claimed result. □ 
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